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Abstract of Thesis Presented to the Graduate School 
of the University of Florida in Partial Fulfillment of the 
Requirements for the Degree of Master of Engineering 
GROUNDWATER MODEL CALIBRATION 
FOR A HYDROCARBON PLUME IN A 
SANDY, SURFICIAL AQUIFER 
By 
STEVEN RAYMOND SCANLAN 
December, 1988 
Chalrman: W. Lamar Miller 


Cochairman: Paul A. Chadik 
Major Department: Environmental Engineering Sciences 


A gasoline-contaminated surficial aquifer, the site of 
four years of remediation efforts, was investigated using a 
semi-quantitative soil extraction procedure to determine the 
vertical distribution of sorbed contaminants present. The 
effectiveness of contaminant flushing in highly contaminated 
soils by aquifer augmentation was demonstrated. The site 
was modeled using BIOPLUME II, an updated version of the 
USGS Solute Transport Model. Although the model was 
effectively calibrated to the groundwater flow through the 
Site, the extreme vertical contaminant distribution and the 
variation of groundwater contaminant concentrations as a 


function of aquifer elevations and time prevented the 


VAIL IL 





acceptable calibration of the pollutant transport model to 


this site. A three dimensional model should be considered. 
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CHAPTER I 
INTRODUCTION 

Groundwater is a valuable resource, and has become 
recognized as such over the past few years. It represents 
96 percent of all fresh water resources in the United 
States, an estimated 36 quadrillion gallons (DeHan, 1981). 
However, its importance, based on use, varies significantly 
across the United States. In 1970, groundwater contributed 
29 percent of the public water supply for the 31 eastern 
states and 46 percent for the 17 western states (Freeze and 
Cherry, 1979). These figures exclude private wells used by 
an estimated 40 million people in the United States, 
particularly in rural areas where most water is consumed 
paWistt..2YGReiicant pretreatment (Miller, 1985). In 
Florida, 92 percent of all drinking water is obtained from 
groundwater sources (Florida Department of Environmental 
Regulation, 1985). 

The awareness of groundwater's importance improved as 
demand for water was increasing, the capacity of known, 
economically available sources was limited, and the quality 
of many known sources was deteriorating from outside 
contamination. The problem with availability of 


uncontaminated groundwater was further exacerbated by 





overpumping of some coastal aquifers, resulting in land 
subsidence and salt water intrusion. 

Potential sources of groundwater contamination are 
numerous, and have become commonplace in our society. They 
include underground storage tanks. Across the country, 
there are an estimated 1,400,000 installed underground 
storage tank systems with capacities greater than 1,100 
gallons, 84 percent of which are made of unprotected steel. 
About 50 percent of these tanks are used for petroleum 
products (Kosowatz, 1988), and the United States 
Environmental Protection Agency estimates that ten percent 
of these petroleum storage tanks are currently leaking. 
Based on benzene being one percent of unleaded petroleum 
products, and the current maximum contaminant level for 
benzene being one part per billion as established by the 
plore. Departient of Environmental Regulation, one gallon 
of gasoline has the potential of contaminating over eight 
million gallons of groundwater (Cherry, 1987). Ina 1978 
study of leaking petroleum storage systems, the Petroleum 
Equipment Institute found 81.6 percent of tank leaks and 98 
percent of piping leaks were the result of corrosion (Kost 
and Parish, 1986). 

Over the past 15 years, significant legislation at all 
levels of government has been enacted to regulate or 
prohibit certain practices commonly linked with groundwater 


contamination. Although there has been progress, 





Significant sources of potential contamination remain, and 
existing contamination plumes continue to spread. 

Once an aquifer is contaminated, the recommended 
cleanup option is dependent upon the characteristics of the 
contaminant, extent of the contaminant plume, physical 
characteristics of the site, and sound engineering practices 
(Nyer, 1985). At best, an aquifer cleanup operation will be 
expensive, difficult, and long term. Considering the large 
investment of time and money for cleanup, predictive models 
of the groundwater system have proven extremely useful in 
determining relative effectiveness of plume control and 
remediation plans (Freeberg et al., 1987). 

There have been three types of predictive models 
commonly applied to groundwater systems: sand tank models, 
analog models, and mathematical models. A difficulty with 
all models is readily duplicating conditions found in the 
field. While the sand tank models and analog models enjoyed 
considerable use prior to the 1960's, they fell out of favor 
when high-speed digital computers became available. While 
mathematical models have been used since the late 1800's, it 
has been computers that led to their widespread use (Wang 
and Anderson, 1982). The advent of the personal computer 
has brought the capability of mainframe computers to the 
desk top. The major limitations of the personal computer, 
in comparison to a mainframe, are limited data set/code 


size, and the speed at which the program is executed (Westly 





and Robertson, 1985). With the recent development of the 
Operating System/2 software, the data set/code size 
restrictions have been reduced, and technological advances 
are increasing the personal computer's operating speed. 
These advances, coupled with the recent proliferation of 
modeling programs for personal computers, have made 
mathematical models virtually essential for all who work 
with groundwater and groundwater contamination problems. 

While mathematical models have simplified the 
predictive process for groundwater contamination assessment, 
transport and remediation, they still require considerable 
Site specific hydrogeologic information to be effective. 
They should also incorporate any biological degradation, 
chemical reactivity, and physical processes such as solute 
adsorption to and desorption from the soil. 

The research described in this thesis concerned a 
gasoline contaminated surficial aquifer at the University of 
Florida, Institute for Food and Agricultural Sciences Citrus 
Research and Education Center (IFAS-CREC), Lake Alfred, 
Florida. Based on previous hydrogeologic characterization 
of the site, the objectives were to: 

1. Identify the vertical extent of adsorbed 
hydrocarbons on the site, 


2. Develop a water balance for the system, 





3. Select a mathematical solute transport model that 
adequately describes the contaminant transport for the field 
site under study, and 

4. Use water table elevations, aquifer 
transmissivities, and groundwater contaminant concentrations 


to validate the selected model. 





CHAPTER II 
LITERATURE REVIEW 
This chapter presents a review of literature pertaining 

to the contamination site with its operational recovery 
system, and the development and application of mathematical 
models. The major areas of discussion are the history of 
the project site, the development of mathematical models, 
the model calibration process, and case histories of model 


application. 


e Projec ite 


The University of Florida, Food and Agricultural 
Sciences, Citrus Research and Education Center (IFAS-CREC) 
is located just north of Lake Alfred, Florida, at 81° 42.9' 
west longitude, and 28° 06.1' north latitude. The facility 
borders on a wetland to the west which eventually drains 
into Lake Swoope, and the site itself is reclaimed wetland. 

The area is underlain by the Hawthorne Formation, a 
clayey-sand confining layer in excess of 11 feet thick. The 
surface of this layer ranges in depth between 6.4 feet and 
12.8 feet beneath the ground surface in the area under 


study. This layer maintains the surficial aquifer in the 





sandy fill material overlying it. The water table surface 
generally follows the ground surface elevation, but varies 
by a few feet depending upon the background flow of the 
aquifer, the existing climatic conditions, and the volume of 
recycled groundwater and fresh water being sprayed on the 
spray field. However, the typical water table elevations 
vary from a maximum of 140 feet above mean sea level on the 
spray field to a minimum of 130 feet near the wetlands. The 
water table is, on average, about four feet below the ground 
surface. Figure 2-1 shows the physical layout of the site, 
and identifies the monitoring and recovery wells. 

Accidental Release and Subsequent Cleanup Efforts 

The following is a summary of Killan's (1987) 
historical overview, readers are referred to his original 
document for a more detailed discussion. 

The University of Florida, Institute of Food and 
Agricultural Sciences, Citrus Research and Education Center 
operated an on-site petroleum refueling system consisting of 
a 7,000 gallon, on-grade gasoline storage tank, a 3,000 
gallon underground gasoline storage tank, a 300 gallon 
underground diesel storage tank, and the necessary 
underground piping and pumps to operate the system. In May 
1983, a loss of leaded gasoline was detected from the on- 
grade storage tank. The tank was removed from service in 
June 1983, and later dismantled for subsequent disposal. 


Air pressure tests of the tank and the underground fuel 
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Figure 2-1. Map of Lake Alfred Study Site. 





transfer lines were performed, indicating leaks, probably 
the result of corrosion. A review of the CREC gasoline 
delivery and pumping records indicated 7,500 to 8,000 
gallons of gasoline were released. 

Between January and May 1984, a series of 12 monitoring 
wells were installed to monitor the contaminant plume and 
facilitate recovery operations. The free product recovery 
operation was conducted between 17 May and 04 June 1984. 
Based on the benzene content of water samples taken 
throughout the operation, an estimated 6,000 gallons of 
gasoline were removed from the aquifer. 

In January 1985, again following a review of their 
delivery and pumping records, IFAS-CREC personnel discovered 
a 50 gallon per day gasoline leak from their 3,000 gallon 
underground storage tank. Subsequently, both underground 
storage tanks were taken out of service in February 1985, 
emptied and filled with sand. 

There has been virtually no free product observed in 
any monitoring well since May i985. 

The drinking water supply, drawn from the underlying 
Floridan Aquifer by a pump housed in building 12, was 
regularly tested for benzene and total hydrocarbons. The 
benzene analysis, performed by a contractor, was 
consistently less than one part per billion. With one 
notable exception, the total hydrocarbon concentration was 


less than one part per million. Upon investigation, the 





10 
exception was found to be the direct result of a quality 
control error in the laboratory. The adjacent wetlands were 
also periodically tested, and a detectable concentration was 
routinely found, but not in sufficient concentrations to 
indicate the contaminant plume ever reached the wetland 
(Rillan, 1987). 

Groundwater Recovery/Spray Aeration System 

The groundwater recovery/spray aeration system depicted 
in Figure 2-1 has been in use since August 1986, except for 
the December 1987 addition of the small intermediate- 
pressure, single nozzle rotating head sprinkler near the 
northwest corner of the electric power substation. Since 
completion of the free product recovery, groundwater 
recovery from select wells, followed by spray aeration has 
continued to strip the volatile organic compounds from the 
water, Agee the less contaminated, oxygenated water to 
the aquifer. 

The operational groundwater restoration system can be 
broken down into three distinct components: the groundwater 
recovery system, the spray aeration system, and the aquifer 
augmentation system. Each component will be discussed 
individually. 

Groundwater recovery system. The groundwater recovery 
system is comprised of submersible pumps (Goulds, Model 
25EL10422) placed in wells RAP-1, UF-2M, and RAP-3. Each 


pump is controlled by a Pump Sentry controller (Model 





13 
AW100), which allows the well to pump dry, shut off for a 


preset time period, and then reinitiate the pump cycle. The 
placement of wells RAP-1 and UF-2M were made to prevent the 
further expansion of the contaminant plume, as well as 
provide recovery points in the areas of highest contaminant 
concentration. RAP-3 serves as a backup well for UF-2M, 
minimizing the threat of release of contaminated groundwater 
to the wetland, as well as providing a back flushing 
capability of the surrounding soil by drawing water from the 
wetland through the surrounding contaminated soil, and 
pumping it into the system. 

For contaminants with a density less than water, the 
use of recovery wells and injection wells as a means of 
hydrodynamic plume control is well documented (Nyer, 1985). 
The effectiveness of these wells is further supported by the 
lack of Pe ie ona plume movement since their operation 
began (Miller, 1988). 

Spray aeration system. The spray aeration system 
consists of 30 low flow rate spray heads installed on 24 
inch vertical risers distributed across the spray field as 
shown in figure 2-1. The droplet size produced is a 
function of the water pressure of the system; at low 
pressures, small droplets are produced, while at high 
pressures, a fine mist is produced. In both cases, the 
air/water interface area is significantly increased for a 


given mass of water. Since the alkylbenzenes, which compose 





12 
most of the gasoline, are highly volatile and have 
relatively high Henry's Law constants (Nyer, 1985), the 
large air/water interface results in a phase change of the 
contaminant from the liquid phase to the gaseous phase. In 
addition, the oxygen concentration gradient at the interface 
will result in oxygen movement from the gaseous phase to the 
liquid phase, thus returning an oxygen enriched, cleaned 
water to the aquifer, when the aerated water infiltrates the 
spray field and percolates to the underlying aquifer. 

Aquifer augmentation system. The aquifer augmentation 
system consists of a double nozzle, medium intermediate- 
pressure sprinkler located in the center of the spray field, 
and the single nozzle, small intermediate-pressure sprinkler 
at the northwest corner of the electric power substation. 
These sprinklers can deliver an additional 15,000 gallons of 
fresh milena Aquifer water per day to the spray field. 

Since gasoline has a density less than water and is 
immiscible in water, the released product will move downward 
under the force of gravity, wetting soil surfaces, with 
minor capillary forces creating some lateral spreading 
during the vertical movement of the plume. Once the 
gasoline reaches the water table, it may create a slight 
depression in the water table surface, but will spread as a 
result of elevation gradients and later, capillary forces, 
with small amounts of hydrocarbons dissolving in the water. 


Free product recovery operations will generally remove 





13 
Significant amounts of the product reaching the water table, 
but some will remain in the soil, either adsorbed to organic 
material (Jury, 1986), or held between soil particles by 
capillary forces. However, since the non-ionic, non-polar 
alkylbenzenes adsorb almost exclusively to organic material, 
and most sand aquifers with a significant water yield 
contain little carbon, Cherry (1987) contends adsorption is 
not a significant process in delaying the contaminant 
movement through such an aquifer. Should the water table 
drop in elevation after the free product recovery is 
completed, the layer of adsorbed and surface tension 
immobilized product will remain in place (Freeze and Cherry, 
1979). These perched contaminated layers are susceptible to 
biodegradation by native bacteria found in most soils, but 
as long as the layer remains, it poses further contamination 
threats to the underlying aquifer because small quantities 
of the contaminant can be released to the aquifer by the 
percolation of water through the contaminated zone, or more 
significant quantities will be released should the water 
table become elevated and engulf the contaminated soil 
region. Since the contaminant is immobilized in the soil, 
the driving force for release of the contaminant is the 
concentration gradient between the contaminated soil and the 
"fresh" water moving through it (Nyer, 1985). While the 
concentration diffusion process is not instantaneous, it is 


a function of the hydrocarbon solubility, the rate of 





14 
transfer from the solid phase to the aqueous phase, and the 
contact time between the two phases. Therefore, elevation 
of the water table or a continuous downward percolation of 
water will have a more significant effect on release of the 
immobilized contaminant than will an infrequent movement of 
water such as that resulting from a small rainfall. 
Additional research on the contaminant transfer from the 
solid to liquid phase under field conditions needs to be 
conducted (Cherry, 1987). 

The aquifer augmentation with the fresh water 
sprinklers elevates the water table, provides a continuous 
flushing action through the soil containing the adsorbed 
hydrocarbons, and increases the groundwater elevation 
gradient across the site. 
Natural Factors Affecting Groundwater Recirculation System 

wtHTe the ideal conditions for modeling a groundwater 
system would be a static situation with neither loss nor 
gain of water, this is an impossible condition in the field 
because of natural occurrences, particularly evaporation, 
evapotranspiration, and precipitation. 

Evaporation. Evaporation is the direct vaporization of 
water from a wet surface. The amount of evaporation 
occurring is a function of relative humidity, air 
temperature, and wind speed. To standardize its 
measurement, the National Weather Service (NWS) has adopted 


standard dimensions, construction materials, and operating 





1) 
procedures for evaporation pans in their weather stations 
(Smajstrla et al., 1982). There is an NWS standard weather 
station operated by IFAS-CREC personnel at the Lake Alfred 
site, and daily measurements of pan evaporation are made. 

Measured pan evaporation values from a properly 
maintained pan, corrected with the appropriate pan 
coefficient, are accurate indicators of the amount of 
evaporation occurring from free water surfaces and wetted 
solid surfaces under similar exposure conditions (Decker, 
19166 ) . 

The use of sprinklers also increases the amount of 
evaporation occurring. The spray emitting from sprinkler 
heads is exposed to the same evaporative conditions during 
its flight through the air, but its increased air/water 
surface area makes it more susceptible to evaporation. 
BRE Lore, the finer the water droplets emitted, the greater 
the evaporation, under the same climatic and physical 
conditions. Water loss due to wind drift must also be 
considered if the drift results in depositing the water 
outside of the study area. 

Experimentation and measurement of evaporation from 
different types of irrigation sprinklers has been conducted 
to reduce evaporative losses. Ali and Barefoot (1981) 
provide an overview of several noteworthy experiments. 
Early pivoting irrigation sprinklers were set with a 


trajectory angle of about 26 degrees above horizontal. 





16 
Although the evaporation was found to be a function of the 
existing climatic conditions, operating pressure, and nozzle 
size, many of these sprinklers were found to have an 
evaporative loss ranging between 15 and 35 percent for a 
7.26 millimeter diameter nozzle to between 40 and 52 percent 
for a 6.32 millimeter nozzle. To reduce the evaporation 
from these sprinklers, their trajectory angle has been 
reduced to seven degrees above horizontal, thus reducing the 
air/water interface time (Ali and Barefoot, 1981). Low 
trajectory irrigation sprinklers are used on the Lake Alfred 
spray field. 

Ali and Barefoot (1981) conducted 72 tests of low 
trajectory irrigation sprinklers under various climatic and 
Operating conditions. Defining evaporative loss as "the 
combined loss of evaporation from the spray and the drift 
loss," they used the following equation for its computation: 


Sprinkler vo e - Collected volume x 100%. 


= Sprinkler volume 


The loss results reported from the tests ranged between 
OQ and 48 percent. The data were analyzed using stepwise and 
least square regression procedures to provide the following 


model: 


Loss = 7.95 - 0.40 RH + 0.83 T + 0.85 WS + 0.03 PR 
+ 25/1) RAT 


where RH = relative humidity (%), 





Wy 


m = temperature ('C), 

WS = wind speed (kilometers/hour), 
PR = pressure (kilopascals), and 
RHT = riser height (meters). 


This model had a square of the multiple correlation 
coefficient (R‘) of 0.91, a standard deviation (S) of 4.42 
percent, and a coefficient of variation of 16 percent (Ali 
and Barefoot, 1981). 

While the evaporative loss of the low pressure spray 
systems are less studied, they are regularly used as 
undertree micro irrigation systems in Florida citrus groves. 
In such applications, the operating pressure is low enough 
to produce small droplets of water. The shielding effect of 
the trees reduce the droplets' exposure to the wind and sun, 
contributing to an average application efficiency of 80 
percent (Smajstrla et al., 1988). 

An additional complication which may occur is the 
evaporative effect of the irrigation system may lower the 
air temperature and raises the relative humidity in the 
immediate area, thus reducing evaporation from the wetted 
ground surface (Decker, 1966). 

Evapotranspiration. Evapotranspiration is the 
combination of evaporation and transpiration, the flow of 
water vapor from the interior of a plant to the atmosphere. 


The amount of evapotranspiration occurring is again 





18 
dependent upon the existing climatic conditions, as well as 
the crop canopy covering the ground. 

While there have been several methods developed for 
predicting evapotranspiration, the easiest and most widely 
used is the pan evaporation method. The equation is 


provided below: 


Evapotranspiration = k xX k, X Pan Evaporation 


C 


where k. = crop coefficient, and 


Ka evaporation pan coefficient. 


Experiments conducted in Fort Lauderdale, Florida 
obtained a crop coefficient of 0.9 for turfgrass. Past 
experiments have also shown that the evaporation pan 
coefficient is variable throughout the year, but under 
Florida conditions, 0.7 is a good year around value (Pair et 
al., 1983 and Jones et al., 1984). 

Precipitation. Precipitation 1s measured at all NWS 
Weather Stations, so this information is readily available. 
Agricultural experts break down precipitation into effective 
rainfall, runoff, and percolation. Effective rainfall is 
the moisture temporarily stored in the soil to meet the 
evapotranspiration requirements of the crops, runoff is 
water lost over the ground surface without ever entering the 
soil, and percolation is the water moving below the crop 


root zone toward the aquifer (Jones et al., 1984). This 
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study will be concerned with effective rainfall and 


percolation. 


Historical Development of Mathematical Modeling 


Mathematical models have been used since the late 
1800's, but 1t was not until the availability of high-speed 
digital computers in the 1960's that their use became 
widespread. 

In 1856, French engineer Henry Darcy published his 
report on the water supply of Dijon, France. Included in 
that report were results of experiments of water flow 
through a sand-filled cylinder. From this, he developed the 
empirical relationship known as Darcy's Law. This is 
commonly attributed as the birth of groundwater hydrology as 
a Ree Cative science (Freeze and Cherry, 1979). 
Groundwater Flow Models 

Groundwater flow models are concerned strictly with the 
flow of water through porous media, the aquifer material. 
They have proven their worth in groundwater planning and 
management. Numerical calculations of Darcy's Law and other 
refinements such as the Theis equation for well drawdown 
have allowed groundwater flow models to become reliable 
predictive tools if calibrated with accurate data for the 
site. The information required includes aquifer parameters 


(porosity, transmissivity, saturated thickness, 
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storativity), physical conditions of the site (initial 
conditions, boundary conditions, locations of pumping and 
injection wells), and a historical database with which to 
calibrate the model. There are several available articles 
and textbooks which discuss these models in detail, 
including Wang and Anderson (1982), and to a lesser extent, 
Freeze and Cherry (1979). 

n t Mode 

By the mid 1970's, groundwater flow models were 
experiencing widespread use and had become largely accepted 
in the groundwater community. This was also the period that 
groundwater contamination problems began receiving 
considerable attention. These factors resulted in research 
efforts being extended to include the modeling of physical, 
chemical and biological processes in groundwater 
contamination. The legal profession often relies on them as 
well, making use of solute transport models in litigation 
cases involving contaminated groundwater (Geraghty, 1984). 

In surface waters, solute transport has been modeled 
using the advective-dispersion equation for many years. 
However, there is a Significant difference between fluid 
flow and solute transport in surface waters and that through 
porous media. While all of the requirements of the 
groundwater flow model still apply, increased dispersion and 
adsorption effects, lower transport velocities, and greater 


aquifer heterogeneities prevent the direct translation of 
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the surface water advective-dispersion equation to 
groundwater transport. Freeze and Cherry (1979, Appendix X) 
provide an overview of the derivation of the advective- 
dispersion equation as it applies to groundwater, while Bear 
and Bachmat (1984) proceed through a rigorous development of 
the basic transport equations. Bedient et al. (1985) 
present an easily understood overview of solute transport 
modeling and the factors which affect it, targeting the 
basic level hydrologist, hydrogeologist, and ground water 
resource manager without a background in modeling. 

Advection. Advection is the movement of solute within 
flowing groundwater. If advection is the sole force moving 
a solute, it will move in the direction of groundwater flow, 
at the seepage velocity of the groundwater, and there will 
be an extremely sharp concentration gradient between the 
leading edge of the plume and uncontaminated water. While 
advection iS a significant force in groundwater, it should 
not be considered the sole source of solute movement 
(Beljin, 1985). 

Hydrodynamic dispersion. Hydrodynamic dispersion (or 
dispersion) is the combination of two known physical 
effects. The first, and most significant effect is 
mechanical dispersion caused by velocity variations in 
magnitude and direction within the pore spaces of the 
aquifer on the microscopic scale, and variations in 


permeability on the macroscopic scale. Molecular diffusion 
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of contaminant from high concentrations to areas of low 
concentration also contributes to dispersion. While 
molecular diffusion will be more significant at lower 
velocities, it is usually considered negligible and is 
disregarded in groundwater flow (Bear, 1979). 


While it is possible to measure in situ dispersion 





coefficients, it is difficult at best, and the results 
obtained have been of questionable reliability. As a 
result, dispersion is normally determined during model 
calibration (Beljin, 1985). 

Biodegradation. It has only been in the past 15 years 
that sizable bacterial populations were detected in the 
ground below the root zone. Prior to this discovery, it was 
assumed to be a sterile region (Wilson et al., 1986). These 
microbial populations have been found to develop in 
mMicrocolonies attached to the soil particle surface. 
Investigators have measured attached microbial populations 
on the order of 10° bacteria per gram of dry soil in 
surficial aquifer material, to a depth of six meters (Wilson 
et al., 1983 and Harvey et al., 1984). 

Harvey et al. (1984) conducted their experiment in a 
sewage-contaminated Cape Cod, Massachusetts aquifer. They 
found no significant difference in the attached microbial 
population in the contaminated plume or in an uncontaminated 
section of the aquifer. However, the population of smaller, 


free-floating bacteria within the contaminant plume had a 





Z3 
strong correlation with the distance from the contaminant 
source. Similarly, in a southern California gasoline spill, 
hydrocarbon bacteria of 50,000 per milliliter of sample or 
more were found in wells containing free product, while a 
nearby uncontaminated well had only 200 bacteria per 
milliliter of sample (Borden and Bedient, 1986). 

While it 1S a growing area of research, there is little 
Known about microbial dynamics and simultaneous 
biodegradation of contaminants in an aquifer system. Under 
favorable conditions, research has shown significant 
reduction of degradable organic contaminants, enough to 
indicate that solute transport models of aquifers under 
these conditions should account for biodegradation. The 
favorable environmental conditions include a degradable 
carbon source, dissolved oxygen, oxidation-reduction 
potential, and the availability of mineral nutrients. Two 
vital requirements for the bacteria are the degradable 
contaminant, serving as a food/energy source, and the 
dissolved oxygen, serving as the primary electron acceptor 
(Widdowson et al., 1987). 

Given favorable conditions and naturally occurring 
aerobic bacteria commonly found in aquifers, benzene, 
toluene, xylenes, and other alkylbenzene components of 
gasoline are readily degraded (Angley, 1987). These same 
compounds have also recently been found to be degraded under 


anaerobic conditions in methanogenic river alluvium 
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contaminated with landfill leachate and a sandy surficial 
aquifer contaminated with aviation gasoline. However, 
anaerobic degradation occurs at a significantly slower rate 
than aerobic degradation (Wilson et al., 1986). 

Aerobic bacterial populations are believed to increase 
until limited by oxygen, degradable organic material, or 
essential nutrients. In most contaminated aquifers, 
dissolved oxygen is assumed to be the limiting factor 
(Wilson et al., 1986). The oxygen consumption in areas of 
high contaminant concentration will be great enough to 
produce significant oxygen concentration gradients. While 
these gradients may induce the diffusion of oxygen from the 
phreatic zone or from oxygenated waters outside the plume 
into the plume (Wheeler et al., 1987), this process is 
unable to keep up with the oxygen consumption rate. 
eeaeuer, there 1S a great deal of uncertainty concerning the 
vertical dispersion coefficient of a conservative tracer. 
Recently measured vertical dispersion coefficients have been 
on the order of molecular diffusion coefficients, orders of 
Magnitude lower than previously believed (Borden and 
Bedient, 1986). 

The rapid oxygen consumption within the plume produces 
an anaerobic zone in the region of reduced contaminant 
concentration, resulting in biomass decay and significant 
reduction, if not total cessation, of aerobic degradation of 


the contaminant. Near the fringes of the plume, the biomass 
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decay appears to be reduced, and this is believed to be the 
result of oxygen diffusion from surrounding, oxygenated 
waters. This oxygen diffusion along the lateral edges of 
the plume is also believed to retard the lateral movement of 
the plume (Widdowson et al., 1987). 

Angley (1987) conducted a series of aromatic 
hydrocarbon degradation experiments using contaminated 
groundwater from the Lake Alfred site. Due to the rapid 
rate of degradation, he concluded the native bacterial 
population was acclimated to the hydrocarbons present. The 
first order rate equation provided the best fit to the 
measured experimental data. 

Card (1987) reviews case histories using biological 
remediation of contaminant plumes, as well as viable methods 
of oxygen augmentation to the oxygen deficient waters of the 
plume. 

Adsorption. Once a contaminant has adsorbed to a soil 
surface, it will come to equilibrium with the surrounding 
fluid, thus gradually reducing the adsorbed contaminant 
concentration. This has the effect of delaying (retarding) 
the passage of the contaminant through the aquifer. For 
compounds undergoing linear instantaneous adsorption, this 
retardation factor (R) may be calculated as shown below 


(Borden and Bedient, 1987): 


mR = Vi40 = 1osepkdi/n 
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where Vv, = ground water velocity of 
nonadsorbing compound, 


V 


A ground water velocity of 
adsorbing compound, 


Kd = partition coefficient (cm'/g 
compound adsorbed), 

P = soil bulk density (g/cm ) and 

n = soil porosity. 


Adsorption may significantly increase the rate of 
biodegradation of a contaminant by allowing oxygenated water 
to overtake the retarded contaminant plume. This is 
particularly true of oxygenated water entering the trailing 
edge of the plume (Wheeler et al., 1987). Without 
Significant adsorption, oxygenated groundwater will be 
unable to overtake a contaminant plume moving at the same 
velocity, nor does longitudinal dispersion appear to have 
much of an effect. In such cases, transverse dispersion has 
the dominant effect (Borden and Bedient, 1986). 

Adsorptive solute transport models may be reduced to 
two general categories, equilibrium models and Kinetic 
models. The equilibrium models assume there is an 
instantaneous equilibrium between the adsorbed solid phase 
and liquid phase solute, while the kinetic models assume 
solutes are simultaneously moving toward and away from the 
adsorbing surface at variable rates and different rate 
constant coefficients (Saleh et al., 1986). The limiting 


process favoring a kinetic model is normally a chemical 
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reaction or a diffusion process. Grove and Stollenwerk 
(1987) provide an overview of groundwater models 
incorporating chemical reactions, as well as a table 
providing an extensive list of models discussed in recent 
literature and their applications. 

While adsorption is assumed to play a significant role 
in the reduction of hydrophobic contaminants in groundwater, 
there is considerable uncertainty about the solute-soil 
interaction process. Two additional processes known to 
occur, producing similar effects, are partitioning and 
absorption. These three processes are often referred to 
collectively as sorption. Although these processes are 
normally considered fully reversible, recent research has 
demonstrated this is not true under all conditions (Angley, 
1957) . 

Aquifer heterogeneities. Varying values of aquifer 
transmissivity will alter groundwater flow lines, thus 
affecting the mechanical dispersion. In surficial aquifers 
of limited thickness, the surface of the underlying 
confining layer may also have a significant effect on the 
groundwater flow patterns. Twenter et al., (1985) found the 
sloping clay confining layer at a spill site altered the 
flow of water at different flow levels of the aquifer. This 
resulted in the displacement of the contaminant plume's 


lateral edge. 





28 

Killan (1987) reported a swale in the confining layer's 
surface between wells RAP-5 and RAP-6. He also indicated 
additional heterogeneous effects from underground utility 
lines, storm sewers, and compression of the aquifer below 
building foundations . The buried steam lines can also 
induce thermal convection currents. 
Further Modelj Developments 

In 1978, the USGS developed their Solute Transport 
Model, a two dimensional, finite difference model for 
nonreactive contaminants. This model incorporated several 
assumptions which are provided below (Konikow and 
Bredehoeft, 1978, p. 4): 

1. Darcy's law is valid and hydraulic-head 

gradients are the only significant driving 

mechanism for fluid flow. 

2. The porosity and hydraulic conductivity of the 

aquifer are constant with time, and porosity is 

uniform in space. 

3. Gradients of fluid density, viscosity, and 

temperature do not affect the velocity 

distribution. 

4. No chemical reactions occur that affect the 

concentration of the solute, the fluid properties, 


or the aquifer properties. 


5. Ionic and molecular diffusion are negligible 
contributors to the total dispersive flux. 


6. Vertical variations in head and concentrations 
are negligible. 


7. The aquifer is homogeneous and isotropic with 
respect to the coefficients of longitudinal and 
transverse dispersivity. 
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The model uses a block centered, uniformly spaced grid 
system with the nodes at the center of the cell representing 
the average value over the area of the cell. 

The transport equation in the model incorporates the 
Method of Characteristics, a method developed to solve 
hyperbolic differential equations, which the transport 
equation closely approximates if advection is the dominant 
mode of transport, a condition existing in most field 
applications. This method has been used successfully in 
numerous field problems without introducing significant 
numerical dispersion. Konikow and Bredehoeft (1978) include 
a detailed description of the program, an explanation of 
each program subroutine, and a printed copy of the computer 
code, written in FORTRAN. 

This model has several advantages. It is a well 
documented and easily understood program that has been 
verified in both laboratory and field tests. There has been 
a menu-driven preprocessor added to the program to make it 
more user friendly, and the program has been updated 
repeatedly, and improvement efforts are continuing (Rifai et 
ol., L967). 

McDonald and Harbaugh (1985) developed a modular, 
three-dimensional finite difference model for groundwater 
flow that is easily modified by the user to fit existing 
site conditions at a new location. The program consists of 


a series of discrete packages, from which the user selects 
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those which apply to the site under consideration, and then 
executes the program. However, this model does not 
incorporate solute transport. 

The USGS Solute Transport Model was modified to 
Simulate the transport and removal of contaminant and 
oxygen, aS well as microbial activities, under the 
assumptions of instantaneous degradation reactions and no 
retardation of the contaminant. Assuming the reaction is 
oxygen limited and there is an instantaneous reaction 
between the contaminant and oxygen, the microbial 
concentration is not necessary, provided it is present in 
sufficient concentration to degrade the contaminant. 
Sensitivity analysis performed on the model indicated 
microbial kinetics did not significantly affect the 
contaminant distribution, thus justifying the instantaneous 
ee one siinpliéication: The development took place at Rice 
University with the support of the United States 
Environmental Protection Agency and was marketed under the 
name BIOPLUME (Wheeler et al., 1987). 

BIOPLUME II incorporates the additional capability of 
including adsorption of the contaminant to aquifer material 
and anaerobic degradation of the contaminant. In both 
BIOPLUME models, there is no time allowance given for 
bacterial acclimatization to the contaminants present; this 
is assumed to have occurred prior to the initiation of 


modeling. There is also an implicit assumption that three 
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molecules of oxygen are required to totally degrade one 
molecule of contaminant (Rifai et al., 1987). 

Mathematical models have also been used to solve the 
"inverse problem," that is, given limited monitoring well 
data, the model is back calculated to produce the likely 
source of the contaminant. This same process may also be 
used to quantify unknown parameters, given sufficient 
information about the contaminant source and monitoring 
wells (Chieh et al., 1985). Yeh (1986) reviews the 
available literature concerning both of these techniques. 

While most of the models discussed concentrate on 
computer codes written in FORTRAN or BASIC, there are 
several proponents of using computer spreadsheets instead of 


these special programs (Olsthoorn, 1985 and Highland, 1987). 


Model Calibration 


The International Ground Water Modeling Center at 
Butler University has established and promulgated a two step 
testing process for groundwater models. The two steps are 
verification and validation. 

Veragicatio 

The verification process has two purposes, to test the 
accuracy of the computer code under known conditions, and to 
ensure the code is fully operational. Both of these 


purposes are accomplished by running the program with a 
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given set of conditions and a known answer obtained either 
from laboratory testing or from an accepted computer model. 
One benefit of using a published program is this should 
already have been completed prior to publication (van der 
Heijde et al., 1985). 

Validation 

Model validation is intended to determine how well the 
model simulates the actual behavior of the system. The 
first method of doing this is known as calibration, and 
involves adjusting the model's physical parameters to obtain 
a good correlation between the model's values and the 
measured field data. The better form of validation is then 
to take the calibrated model, and use it in a predictive 
mode beyond the limits of the time period of the data set 
used for calibration (van der Heijde et al., 1985). 

This calibration process is normally done manually, 
using trial and error. However, there have recently been 
some numerical estimating processes developed for parameter 
estimation (Yeh, 1986). One key point concerning 
groundwater models is that there is not one unique set of 
system parameters for the calibrated model. This is clearly 
demonstrated in Freyberg's (1988) article. Another 
interesting point brought out in the article is the best fit 
model was the one utilizing a few large regions in the 


conductivity field, while the worst fit model used numerous 
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small regions in the conductivity field, attempting to 


"tweak" the model to match the data. 


Case Histories of USGS Model Application 


At the Coast Guard's request, the U.S. Geological 
Survey investigated an organic solvent contamination plume 
in a surficial aquifer in East Bay Township, Michigan. The 
long narrow contaminant plume discovered indicated oxygen 
diffusion into the plume along its lateral edges had 
significantly retarded the lateral spread of the plume. The 
model was calibrated using measured field data, and was then 
used as a design tool to develop the most effective recovery 
well system to prevent further movement of contaminant from 
the site (Twenter et al., 1985). 

The USGS model was also applied to a creosote waste 
site in Conroe, Texas. On the site, there is a hydrocarbon 
plume contained within a larger chloride plume. The 
chloride is believed to have originated as a contaminant or 
decomposition product of pentachlorophenol, a compound used 
in wood treatment processes and present in the waste lagoon. 
The chloride concentrations and distributions were used to 
calibrate the model initially, and then the model was 
applied to the hydrocarbon plume. The original computer 
code had to be modified to incorporate the biodegradation. 


These modifications were the basis of the BIOPLUME program. 
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They provided a good fit of the field measured values 
(Borden et al., 1986). 

Freeberg et al. (1987) applied the United States 
Geological Survey (USGS) Solute Transport Model to a 
trichloroethylene-contaminated, shallow, sandy aquifer. Due 
to curved flow lines on the site's potentiometric surface, 
the configuration and values for a constant head boundary 
were adjusted during model calibration. They also found 
varying transmissivity spatially across the grid produced 
better results, and the transmissivity values used in the 
model varied by 23 percent from the measured value. The 
model was calibrated by minimizing the error (E) between the 
predicted and measured contaminant plumes, using the 


following formula: 


Ese |P, - | 


where rE = predicted concentration of 
TCE at well x, and 
M, = measured concentration of TCE 


at well x. 


In an experiment to determine in situ retardation 





factors for naphthalene and paradichlorobenzene, Borden and 
Bedient (1987) used three wells in a line, six feet on 
center, with the center well used as an injection well, and 
the two outer wells were pumped to induce a greater 
hydraulic gradient and monitoring wells. They attempted to 


calibrate this system using the USGS Solute Transport Model, 
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however were unable to obtain an acceptable solution with 
either uniform or non-uniform permeability distributions. 
This case illustrates the reason that the code authors do 
not encourage its use in problems dominated by radially 
convergent and divergent flow patterns around wells. 
However, El-Kadi (1988) evaluated three situations of such 
conditions, finding two produced acceptable results by 
utilizing a finer mesh. Efforts are presently underway to 


modify the code to improve this aspect of the model. 





CHAPTER III 
MATERIALS AND METHODS 
The efforts of this study can be divided into three 
Main categories: semi-quantitative soil core analysis, water 
balance of the system, and validation of a solute transport 


model at the Lake Alfred site. 
Soi € lysis 


To gain an appreciation for the vertical distribution 
of sorbed hydrocarbons in the aquifer, soil cores were taken 
from several locations on the site. These cores were taken 
by three different methods over the course of this study. 
Coring Procedure 

The initial method of coring employed a trailer-mounted 
hydraulic sampler. Although the equipment functioned 
properly, the cores were difficult to retrieve from the core 
barrel, and impossible to retrieve intact. A split spoon 
sampler would have been preferred. In addition, reference 
marks were not used to differentiate between the undisturbed 
Sample and collapse material from the bore walls following 
removal of the previous core, which became particularly 


Significant once the aquifer was entered. These problems, 
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along with questions about the integrity of the samples 
resulted in the data being discarded. 

The second coring method used a stainless steel, 
manually operated bucket auger. While most problems 
encountered with the first method were solved, the continual 
collapse of surrounding aquifer material into the bore hole 
following the removal of the auger prevented this from being 
a viable option, particularly since the upper portion of the 
elevated aquifer was of greatest interest. 

In the interest of sample integrity, the final method 
was developed and used. Since the area of interest was the 
upper portion of the aquifer, the initial two to three feet 
of overburden were excavated with the manual bucket auger 
and discarded. An eight foot length of 14 inch (14 inch 
inner diameter ) LEXAN tubing, sharpened at one end, triple 
rinsed with 2-propanol and water, was then placed into the 
existing hole, the length remaining above ground was 
recorded, and it was driven to the desired depth by using a 
fabricated six-pound sliding well-point hammer. A steel 
plunger with a rubber gasket was inserted into the open end 
of the tube, gasket first, and pushed down until it 
contacted the sample. A large diameter dowel rod was 
inserted to occupy the void, and the end of the tube was 
Plugged with a number 6 rubber stopper. The plunger, dowel 
rod, and stopper minimized the loss of core material due to 


the suction effect of the aquifer upon tube removal. The 
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core was then driven from the tube by a horizontal force 
through the open end, or through controlled, sharp vertical 
movements, thus allowing the mass of the plunger to drive 
the core out. The coring assembly is shown in Figure 3-1. 

The cores were divided into discrete samples of 
measured length, with each sample placed in a sealed "zip 
lock" sample bag, marked, sealed in a second "zip lock" bag, 
and stored on ice for transport to the laboratory. 

actio 

Within ten hours of taking the first core, this 
extraction procedure was performed. While still in the 
sealed bags, the samples were thoroughly kneaded to produce 
a homogenous mixture. A quantity of the soil mixture, about 
15 grams, was placed in a clean 40 milliliter VOC vial, 
weighed, filled with pesticide grade n-hexane, sealed with a 
TOC GH lined septum cap and shaken by wrist action for 10 
seconds. Although it would have been preferable to extract 
the hydrocarbons without headspace in the sample vials, this 
was impossible because of the air entrained in the soil 
Matrix. Instead, the final head space was minimized through 
careful handling and sealing of the vials. The vial was 
again weighed to obtain the mass of hexane added. These 
prepared samples were then stored at room temperature for 
later analysis. 

Based upon the initial site investigation conducted by 


Killan (1987), the contaminant plume was found to contain 14 
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aromatic hydrocarbons in sufficient concentrations to be 
readily identified through gas chromatography. Since the 
groundwater was tested monthly for these analytes, they were 
also used to analyze the soil extracts. The 14 analytes 
were benzene, toluene, ethylbenzene, meta-xylene, 
para-xylene, ortho-xylene, isopropylbenzene, propylbenzene, 
3-ethyltoluene, 4-ethyltoluene, 1,3,5-trimethylbenzene, 
2-ethyltoluene, 1,2,4-trimethylbenzene, and 
1,2,3-trimethylbenzene. 

Aromatic standard preparation. The stock standard was 
prepared by adding a known mass of each of the above 
analytes to pesticide grade methanol to produce a final 
volume of ten milliliters. Because of co-elution of 
analytes, about 0.010 grams of each of meta-xylene, para- 
xylene, 3-ethyltoluene, and 4-ethyltoluene were added, while 
about 0.020 grams of each of the remaining ten analytes were 
used. To produce a calibration curve, the stock standard 
was diluted in pesticide grade n-hexane using dilution 
Boeroseor 1:2,000, 171,000, 1:250;7°12125, 1:62.5, 1740, and 
1:20. All prepared standards were stored in sealed two 
milliliter serum vials in an explosion-proof freezer until 
required. Storage of these standards did not exceed 30 
days. 

Soil extract analysis. The soil extracts were analyzed 
on a Perkin-Elmer 8410 Gas Chromatograph having a 30 meter 


DB-1 capillary column with a 0.53 millimeter inner diameter 
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and a three micrometer coating. All samples were injected 
directly through a Thermogreen LB-2 11 millimeter septum, 
with the septum being replaced after every ten injections. 
Each sample was analyzed at least twice, and two diluted 
standards were run after every ten sample runs. All extract 
samples were taken by syringe through the septum. Once the 
septum was pierced, the analysis of that sample was 
completed as soon as this method would allow. The 
instrument's analytical conditions are provided in Table 
3-1. 

Recovery study. To verify extraction efficiency of the 
hexane extraction method for the aromatic compounds of 
interest, two recovery studies were conducted on 
representative soil samples obtained from the site. The 
first recovery study was conducted at low concentrations 
without ee head space, while the second study used 
higher hydrocarbon concentrations and included head space 
analysis. In both studies, a Known mass of soil was placed 
in a clean 40 milliliter vial, injected with a measured 
volume of stock standard, sealed with a Teflon lined septum 
and cap, and allowed to equilibrate at room temperature for 
20 minutes. Equilibrium was assumed to have been reached 
within this 20 minute period. In the first study, the 
sample vials were opened, filled with pesticide grade 
n-hexane, resealed, shaken, and later analyzed by the 


procedures discussed previously. In the second study, 
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TABLE 3-1. Perkin-Elmer 8410 Analysis Conditions for 
Soil Extracts. 


Oven Temperature (°C) 35 10 280 
Iso Time (minutes) ooo 15%0 20 
Ramp Rate (°C/minute) 5.0 30.0 

FID Sensitivity High 

Detector zero On 

Injector Temperature 2500¢C 

Detector Temperature 300) °c 


Flow Rate of Helium Carrier Gas 10 milliliters/minute 
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following the equilibration period, a ten microliter sample 
of the vial's head space was taken and analyzed using the 
aforementioned procedure. The head space analysis was 
repeated before the vial was opened and hexane added. 

The results of this recovery study are included in 


Chapter IV. 


System Water Balance 


To obtain an estimate of evaporation from the spray 
heads, an experiment measuring the amount of water reaching 
the ground surface in a known area under two application 
pressures was made. The equipment was setup as indicated in 
Figure 3-2. The pumped volume for each application was 
metered, and the spray heads were identical full-circle 
heads. The collected water was drained from the white vinyl 
collection sheet and stored in small-mouth plastic 
containers, to minimize evaporation, for later measurement. 
Upon completion of each test, approximately 20 minutes, the 
volume of collected water wasS measured using a graduated 
cylinder. 

Daily pan evaporation and precipitation data were 
obtained from the Lake Alfred weather station. Using this 
information, the records of the aquifer augmentation and 
recovery system, and the sprinkler efficiency rates, an 


approximate water balance was calculated. 
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Model Validation 


While the following two sections are presented as two 
discrete topics, the inclusion of the Method of 
Characteristics in the BIOPLUME II code necessitates the 
particle movement calculations be completed, even if there 
is no contaminant present. Therefore, to minimize the 
required computer time, both processes were conducted 
Simultaneously, with the initial emphasis being placed on 
hydraulic flow. Once this was calibrated, the emphasis was 
placed upon the solute transport calibration. 

Groundwater Flow 

A water table database was developed through periodic 
measurements of the water elevation in each monitoring well. 
This information is presented in Appendix A. Metered 
volumes ae water were injected into select wells, to 
evaluate the effect on water table elevation and recovered 
volume of groundwater. The injection and recovery 
information is presented in Appendix B. Based on these data 
and measured aquifer parameters (Appendix C) from Killan 
(1987), a trial and error method was used to calibrate the 
groundwater flow portion of the BIOPLUME II model. 

Three periods of stable system operation were selected 
for the calibration procedure, each representing a different 
aspect of the system: steady flow through the system without 


injection, injection at the upgradient edge of the study 
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area in RAP-5 and RAP-7, and injection in the middle of the 
study area through P-6 and P-7. The model's steady state 
mode was employed, using one week pumping periods and the 
average metered injection and pumped volumes over that 
period. 

The error minimization was performed using this 


Freeberg et al. (1987) equation as a guideline: 


B=. (Be = M, | 


where Er = predicted value at well x, and 


at 


measured value at well x. 


The water elevations in eight wells (M-1, OHM-3, OHM-4, P-5, 
RAP-4, RAP-12, RAP-13, and UF-1E) were used. 

A K-V Associates' Groundwater Flowmeter (Model 30) was 
used to obtain the groundwater flow velocity and direction 
in five selected two-inch monitoring wells. The field 
procedures used were those specified in the system's 
operations and maintenance manual (K-V Associates, Inc., 
1982). To obtain values under static conditions, the 
measurements were taken after the recovery wells and spray 
field had been turned off for over 30 hours. The flow meter 


had been calibrated in a sandy soil matrix prior to field 


use. 
Solute Transport 


The groundwater on site was sampled monthly, in 


compliance with the contract between IFAS and the Department 
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of Environmental Engineering Sciences. The hydrocarbon 
analysis was performed on the Perkin-Elmer model 8410 gas 
chromatograph utilized for the soil extract analysis. In 
addition, a Tekmar LSC/ALS purge and trap system was used to 
concentrate the water samples prior to analysis. The 
program settings are provided in Table 3-2 (Angley, 1987). 

Once the hydraulic flow was calibrated, those 
parameters remained constant, and the solute transport 
Parameters were adjusted by trial and error. The Freeberg 
et al. (1987) equation was applied to the meta-xylene and 
para-xylene concentrations in wells OHM-3, OHM-4, P-5, 


UF-2M, and UF-3W. 


TABLE 3-2. Perkin-Elmer 8410 Analysis Conditions for Water 


Samples. 

1 Z 3 4 
Oven Temperature (°C) 50 70 94 200 
Iso Time (minutes) 520 720 0.0 0.0 
Ramp Rate (°C/minute) 320 3.0 30.0 
FID Sensitivity High 
Detector Zero On 
Detector Temperature 300 °c 


Flow Rate of Helium Carrier Gas 5.5 milliliters/minute 





CHAPTER IV 
RESULTS AND DISCUSSION 


The information in this chapter is divided into four 
sections; the first (The Site) includes field work and soil 
extraction results, the second (System Water Balance) 
presents the initial assumptions, summarizes the 
calculations and discusses the results, the third (Model 
Selection) discusses the model selection criteria, and the 
fourth (Model Validation) discusses the process used and 


problems encountered with the model calibration. 


The Site 


A number of problems arose during the continuous 
operation of the recovery system. 
Operation and Maintenance of Recovery System 

Because of the geographic separation of the project 
Site from the University campus, it was not possible to 
maintain a continuous on-site presence. IFAS-CREC support 
was enlisted to record daily readings, perform emergency 
repairs, and ensure proper day-to-day operation of the 


system. 
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One significant problem encountered was fouling of the 
pump intake screens by bacterial growth in RAP-1 and UF-2M. 
This phenomenon did not occur in RAP-3. The pumps were 
routinely removed from the well and cleaned. After 
returning to operation, the pump would experience flow 
restriction problems after about two months, and would plug 
after three months. This was overcome by weekly dosing the 
well with one gallon of five percent sodium hypochlorite 
solution. With the pump turned off, the hypochlorite 
solution was added to the water in the casing and allowed to 
contact the pump and well casing for about 20 minutes before 
the pump was returned to operation. This periodic dosing 
kept the biofouling problem under control. 

Each pumping well and the irrigation sprinklers in the 
spray field had a dedicated, in-line water meter measuring 
the flow volume. These meters were standard multi-jet flow 
meters, With plastic measuring cylinders. While these 
meters were not new when installed in the system, they 
experienced a 100 percent failure rate between the 
eighteenth and twenty-fourth month of system operation. 

Finally, the spray heads in the spray field were not 
operating efficiently at low-flow periods. This was 
overcome through reduction in number of operating spray 
heads, and reducing some full circle spray heads to one-half 
or one-quarter circle heads. These measures effectively 


increased the operating pressure of the system, producing a 
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finer spray, thus increasing the air-stripping ability of 
the system. 

Groundwater Flowmeter 

A graphic summary of the results from the field 
measurement of the groundwater velocity and direction is 
provided in Figure 4-1, with the numerical values listed in 
Table 4-1. 

Killan (1987) speculated on the aquifer compression 
experienced beneath building 10, producing an effective 
barrier to groundwater flow, diverting it to the south for 
passage between buildings 10 and 12. The test conducted in 
well M-3 supports this idea with its northwesterly 
direction. There is a similar effect, although not as 
dominant, in well OHM-4. The building induced compression 
and the natural movement of water toward RAP-1, even without 
pumping, helps explain why the gasoline pooled south of the 
wash rack, without spreading north to RAP-2 or west of 
buriding 31. 

The pump house (building 12) has a 79 inch deep 
concrete wet-well in its south east corner. The physical 
presence of it would have prevented the passage of free 
product and severely restricted the flow of dissolved 
organics through this area, since the confining layer 
surface is only eight feet below grade at this location. 
The flows in RAP-4 and RAP-6 toward the opening between 


buildings 10 and 12, and the flow of RAP-11 away from this 
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Table 4-1. Groundwater Flowmeter Results for Select Wells 


Computed Velocity Direction? 
Well Reading (feet/day) of Flow 
M-3 5 es 13 298° 
OHM-4 90.6 2.6 2575 
RAP-4 259 >4.9? 305° 
RAP-6 159 4.8 293° 
RAP-11 6.7 <0.9° 149° 


1 360° equals North 


‘ The instrument's calibration curve upper limit for these 
tests was 4.9 feet per day. 


* The instrument's calibration curve lower limit for these 
tests was 0.9 feet per day. 
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opening indicate the existence of a groundwater divide 
between these two areas, probably induced by a combination 
of the pump house and variations in the confining layer 
surface. 

During the installation of additional monitoring wells 
for further research, a shallow swale was discovered between 
buildings 10 and 12, running toward RAP-1. 

Soil Core Analysis 

The purpose of this portion of the study was to 
qualitatively determine the vertical distribution of sorbed 
contaminants in the soil. While the first two coring 
methods proved unsuccessful, the final method produced 
acceptable results. 

Coring process. The LEXAN tubing was an effective 
coring device, withstanding the forces of driving without 
Significant damage. However, it should be noted that the 
site material was sandy-clayey soil. The maximum core 
length that could reasonably be extracted from the core 
barrel by a horizontal driving force was 24 inches, but a 
30 inch length of aquifer material was easily recovered by 
controlled, sharp vertical movements of the core assembly. 

Analytical process. Because of the variation in 
response of the flame ionization detector to each analyte, a 
separate calibration curve was developed for each of the 12 
eluting peaks of interest. There was a minor interference 


between the benzene peak and a small peak attributed to an 
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impurity in the hexane. Since the volume of Sample injected 
was kept fairly constant, about 3.0 microliters, this was 
accounted for in the calibration curve. 

The range of each calibration curve varied between 5 
and 350 nanograms of analyte. The curves may be found in 
Appendix D. From these curves, the mass of analyte was 
determined, given the integrated response of the 
chromatogram. The mass of sorbed analyte in the soil sample 
was calculated from the following equation: 

Meorted = ats Vol, pan, 


where M, = mass of sorbed analyte, 
orbed : 
nanograms solute/gram soil, 


M, = mass of analyte, nanograms, 

M, = mass of hexane in vial, 
grams, 

M, = mass of soil, grams, 

Vol, = volume of injected sample, 


microliters, and 


density of hexane, at 20°C 


Den, 
0.6603 grams/milliliter. 


The method and instruments used were capable of 
detecting lower concentrations than the low limits adopted 
in this study. The calibration curves of several analytes 
began displaying a nonlinear relationship at levels below 
these limits. The high limit of the curve was set so as to 


avoid overloading the column. The curves utilized were 
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therefore restricted to the region of linear relationship, 
without threatening column contamination. If the analyte 
was not detected by the analytical method, it was reported 
as ND and given a sorbed mass of zero, while if it was 
detected but in insufficient concentration for the limits of 
the calibration curves, it was reported as BDL and given the 
sorbed mass value produced by one-half of the lower limit 
value of that analyte. 

Discussion of results. The sorbed hydrocarbon 
concentrations as a function of elevation are plotted for 
each core taken, as Figures 4-2 through 4-6. The samples 
ranged in length from six to nine inches, and the midpoint 
elevation of each discrete sample was used for these plots. 
The tabulated raw data may be found in Appendix E. 

Although the less contaminated cores (Figures 4-3, 4-4, 
and 4-5) appear to vary without a dominate concentration 
profile, this is attributed to analytical results at or 
below the calibration curves' lower limit, and slight 
differences in the soil mass used in the extraction. The 
highly contaminated soil profiles (Figures 4-2 and 4-6) 
indicate a narrow, horizontal region of contamination. The 
contaminated region in Figure 4-6 was reported by Killan 
(1987) to have contained floating free product prior to the 
recovery operations. 

While the lower sorbed concentrations found at the 


upper and lower edges of this region may indicate the result 
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of thoroughly mixing the soil from each discrete sample 
prior to extraction, they may also reflect contaminant 
volatilization, transport within the capillary fringe, 
flushing effects of water movement through the region, or 
biodegradation. More likely, it may be a combination of any 
or all of the above processes. 

Since the free product was known to have pooled east of 
building 31 and south of the wash rack, this is a likely 
area of significant sorbed contaminant. Being covered by 
asphalt, contaminant volatilization and oxygen diffusion 
from the atmosphere into the soil is prevented. 

Angley (1987) reported organic carbon constituted 0.015 
percent of the Lake Alfred aquifer material. In laboratory 
leaking column and equilibrium batch isotherm experiments, 
he measured this aquifer material's retardation factor 
ranging from 1.66 to 1.99 for both meta-xylene and para- 
xylene. Based on these results, a low retardation factor 
(about 1.8) would be expected in this aquifer under 
equilibrium conditions. The presence of the highly 
contaminated soil located just above the saturated zone, and 
the adjacent highly contaminated soil just below the upper 
limit of the saturated zone in the elevated aquifer, should 
serve aS a continuous contaminant source until full 
equilibrium with the surrounding aquifer is reached. The 
water flowing through the contaminated area rapidly becomes 


saturated with hydrocarbons, thus preventing the transport 
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of the total sorbed mass at the seepage velocity of the 
groundwater, which would indicate a retardation factor of 
1.0. The situation is further complicated by small vertical 
dispersion coefficients. This is the type of problem which 
has encouraged the development of three-dimensional solute 
transport models. 

All the monitoring wells on site have screen sections 
fully penetrating the aquifer, and the wells are purged of 
three casing volumes of water before a sample is taken for 
hydrocarbon analysis. This produces a composite sample over 
the depth of the aquifer. 

Recovery Study 

The first recovery study was unsuccessful. This 
failure can be attributed to several factors. First, the 
soil sample used for both studies was slightly contaminated 
with sorbed hydrocarbons. Although the control sample was 
analyzed twice for each study, several detected analytes 
were at or below the lower limit of the calibration curves. 
In such a case, the mass of solute on the column was 
estimated by using one-half the calibration curve's lower 
limit, thus introducing additional errors. Low 
concentrations of the analyte spike and extensive 
volatilization of hydrocarbons with the addition of hexane 
(visible vapors were seen escaping), further complicated the 
procedure. The results are tabulated in Table 4-2. While 


the values in this table may appear reasonable, the 





Table 4-2. Results of Initial Hydrocarbon Recovery Study 


Soil Mass of Total 
Sample Repetition Hydrocarbons Percent 
Number —__Number _ Injected (mq) Recovered 
il a On127 48.3 
2 Onr2 7 3159 
3 OvirZ7 28.4 
a dl 0.637 66.8 
2 0.637 Jom 
3 0.637 95.6 
1 al 1.274 64.1 
2 1.274 67.4 
3 1327 4 81.8 
Z al O2i27 14.5 
2 O.127 -39.0 
3 OviZ7 -28.5 
Z al 0.637 Died 
Z On 6m 59:3 
3 0.637 48.2 
2 1 1.274 Hi 6 
2 1.274 1a 


3 Aiea 13.2 
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variation between replicate analyses of the same sample 
becomes more significant, particularly with regard to 
specific analytes. These tabulated raw data may be found in 
Appendix F. 

The second recovery study included higher analyte spike 
concentrations, and analyzing the sample vial head space 
prior to the addition of hexane. The mass of hydrocarbons 
lost due to the hexane displaced head space was calculated 


uSing the following equation: 


Mead = MxM x (1000 ul/ml) 


Vol; x Den, 
where M,,, = mass of volatilized analyte, 
nanograms solute, 
M, = mass of analyte, nanograms, 
M, = mass of hexane in vial, grams, 
Vol; = volume of injected sample, 


microliters, and 


density of hexane, at 20°C 


Den, 
0.6603 grams/milliliter. 


The mass of hydrocarbons displaced from the vial, when 
combined with the total sorbed hydrocarbons in the sample, 
adequately accounted for the mass of hydrocarbons in the 
spike and that originally in the sample. The results of 
this study are presented in Table 4-3. 

While sufficient repetitions were not conducted to 


establish a rigorous quantitative method, the results are 








Table 4-3. 
Mass of Total 
Hydrocarbons 


Injected (mg) 
2.517 


R259 


a 39 


Repetition 


—Number__ 


1 


2 


Extraction 
Percent 


Recovered | 


695 
O20 
66.7 
62.0 
Sh) 
TAG 
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Results of Second Hydrocarbon Recovery Study 


Total 
Percent 


Recovered 


112 
104 
ogo 1 
92.4 
$2.0 
76 
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sufficient to support the semi-quantitative application used 
in this study. 

Fresh Water Injection. 

During this study, metered volumes of water were 
injected into selected wells to evaluate this method's 
effectiveness in elevating the water table. The injection 
wells, and their period of use are presented in Table 4-4. 
The daily volumes of water injected into each well were 
normalized to a 24 hour period between 0730 and 0730 the 
following day. These normalized volumes are included in 
Appendix B. 

Elevations of the water table surface in selected wells 
as a function of time are presented as Figures 4-7 and 4-8. 
Both figures demonstrate the hydrodynamic influence high 
volume injection wells have on the aquifer, particularly 
between days 727 and 758 (injection in RAP-5 and RAP-7), and 
between days 814 and 832 (injection in P-6 and P-7). They 
also provide an indication of the extent of the regularly 
occurring fluctuations. 

Figure 4-9 presents the combined daily volume pumped 
from the three recovery wells. The general shape of this 
curve is similar to the water table elevation curves, with 
the greatest similarity exhibited by wells closest to the 
recovery wells. It should be noted this is a six day moving 
average curve, thus introducing a delay in exhibiting the 


effects shown in the water table elevations. The low points 





Table 4-4. 


Well 


Injection 


30 
30 
30 
26 


26 


Well Utilization 


Dates Used 
JAN - 24 APR 
JAN - 24 APR 
JAN - 29 APR 
APR - 14 MAY 
APR - 06 MAY 


1988 
1988 
1988 
1988 


1988 


727 
12) 
20 
814 


814 


—eSt2Z 
- 812 
= ey 
ao s2 


- 824 
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Project Days 
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71 
at project days 733 and 813 reflect power outages to UF-2M; 
day 793 reflects the flow meter failure at UF-2M; and days 
750 and 761 reflect disruption of the RAP-1 distribution 
system by an on-site contractor. Similarly, the low points 
at day 841 coincide with the cessation of injection at P-6 
and P-7, and the irrigation sprinklers being secured for two 
days. The final low point at day 854 was caused by the 
irrigation sprinklers being secured for a 12 day period. 

The flushing capabilities of injection wells were also 
demonstrated in this study. The effect was most pronounced 
in the highly contaminated area around P-5, OHM-3, OHM-4, 
and UF-2M. Figures 4-10 through 4-13 present the combined 
meta-xylene and para-xylene concentrations for each of the 
wells. Background data from August 1987 to January 1988 is 
also included to emphasize the effect of the injection 
wells. | 

In P-5, a relatively uncontaminated well in comparison 
with the other three wells, there was a marked decrease in 
contaminant concentration following injection in P-6 and 
P-7. While this concentration has remained suppressed for 
three monitoring periods since the termination of injection, 
continued monitoring is necessary to evaluate the extent of 
recontamination from up gradient areas. 

In OHM-3 and OHM-4, the concentration increased with 
the injection operation in RAP-5 and RAP-7, but considering 


the variation in concentrations over the preceding months, 
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and the decrease in UF-2M's concentration, it is not 
possible to attribute this increase strictly to the 
injection of water. However, following the injection in P-6 
and P-7, the three wells (OHM-3, OHM-4, and UF-2M) recorded 
their highest concentrations in the 2.5 year history of the 
Environmental Engineering Sciences monitoring program. This 
supports the hypothesis that high levels of contaminant are 
present under the asphalt surface which were previously not 


being flushed by the unaugmented groundwater flow. 


System Water Balance 


Sprinkler Evaporation Test 

The sprinklers were not protected from the wind during 
the sprinkler evaporation experiment, in order to simulate 
their normal operating conditions. While the wind was not 
blowing during the 0000 or 0600 tests, it was blowing 
lightly during the 1800 test, and was a significant factor 
in the 1200 test. The results are presented in Table 4-5. 

Low pressure. In the low pressure test, the spray head 
was producing a continuous 12 foot diameter spray pattern 
and the collection sheet covered a 61.3 degree arc, 
extending the full spray radius of the sprinkler. The spray 
head was also isolated from possible overspray from 
surrounding heads. The collection percentage was greatest 


early in the morning, closely followed by the midnight 





Table 4-5. Sprinkler Evaporation Collection Results 


Low Pressure 


Volume Volume Percent 
Time Pumped al Collected (gal) Collected 
1800 88.0 2.88 So All 
0000 96.4 3.93 4.08 
0600 89.3 35 8S 4.32 
1200 96.5 1.29 peace 

High Pressure 
Volume Volume Percent 


Time Pumped (gal) Collected (gal) Collected 


1800 IVS .2 1.49 120 
0000 Ma tale 4 ies 1.24 
0600 OE O pS 30 


1200 Ore al cial ic 
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collection period. While the smallest percent collected 
occurred during the noon sample, as expected, this was also 
the period the wind had its greatest effect on the spray 
pattern, and solar radiation and air temperature were 
nearing their daily peaks. Therefore, it is impossible to 
evaluate from these limited data the extent of evaporation, 


and the extent of loss due to the wind affected distribution 


patterns. 
High pressure. The high pressure test produced an 


intermittent 24 foot diameter spray pattern and was not 
protected from overspray. This collection sheet only 
covered the inner 6.5 feet of a 64.2 degree arc through the 
spray pattern. The results from this phase of the 
experiment show little deviation in the percent collected at 
the designated times. This apparent inconsistency is 
attributed to the unequal spray distribution across the 
pattern's radius, and overspray from surrounding heads. 
Spray System Efficiency Rates 

The low pressure system consists of 15 spray heads on 
the north half of the spray field. Since the prevailing 
wind is from the north, any wind "loss" will still fall on 
the spray field. Unlike the watering systems used in citrus 
groves where the average application efficiency is 80 
percent (Smajstrla et al., 1988), this spray is fully 
exposed to the wind and direct solar radiation. Therefore, 


an application efficiency of 70 percent was applied to the 
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volume pumped from RAP-1 and distributed on this section of 
the field. 

The main high pressure system distributes the flow of 
UF-2M through ten spray heads on the south half of the spray 
field. The heads produce a mist, readily affected by the 
prevailing wind which transports considerable quantities to 
the asphalt road adjacent to the south edge of the spray 
field, resulting in a loss to the system. The fine spray 
and wind loss combine to produce a lower than average 
application efficiency. Therefore, a 50 percent efficiency 
was applied to the flow from UF-2M. 

The intermittent flow from RAP-3 is applied to the 15 
spray heads of the low pressure system. During this period, 
this section also becomes a high pressure system. However, 
because of its location on the spray field's northern half, 
it is not subjected to the high wind loss out of the system. 
Therefore, the flow from RAP-3 was assigned an application 
efficiency of 60 percent. 

The fresh water irrigation sprinklers in the spray 
field are operated at a relatively low pressure, thus 
producing a steady stream and large droplets of water. 

Since Ali and Barefoot (1981) reported evaporative loss 
ranged from 0 to 50 percent, an efficiency of 80 percent was 


applied to this flow stream. 
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Addition Water Sources/Si 

Evapotranspiration. Under identical climatic 
conditions, a wetted surface is assumed to experience an 
evaporation rate equal to the pan evaporation rate 
multiplied by its pan coefficient. In the spray field, the 
constant spray will increase the relative humidity and 
decrease the air temperature, both factors reducing the 
evapotranspiration rate. However, considering the limited 
extent of the spray field, 5,250 square feet, and the 
prevailing wind blowing during the peak evaporation periods, 
it was assumed there is no reduction in the 
evapotranspiration rate of the field. 

Turfgrass in Florida will transpire an annual average 
of 63 percent of the measured pan evaporation (Jones et al., 
1984). The grass covered surface within the area of 
interest is approximately 14,200 square feet, excluding the 
wetted area of the spray field. 

Injection wells. During this study, fresh water was 
injected into the aquifer through select monitoring wells. 
It was assumed there were no evaporative losses from this 
augmentation system. In the unique case of RAP-7, injecting 
adjacent to the model-designated no flow boundary, it was 
assumed 50 percent of this freshly injected water crossed 
the boundary and was lost from the system 

Precipitation. The site has roof gutter systems and 


storm sewer systems which are capable of effectively 
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removing precipitation from improved areas before it can 
enter the soil. The assumptions listed below were made 
concerning this potential water source. 

1. Surface run-off of rain falling on the grassed 
areaS was assumed to be 10 percent of the total 


volume, and 


2. There was no runoff from the developed areas 
onto the grassed areas. 


Groundwater flow. The naturally occurring groundwater 
flow through the area should also be considered. Killan 
(1987) estimated this flow to be 5,500 gallons per day. 
During this study period, a relative decrease in groundwater 
flow was indicated by the low levels of water in the 
adjacent wetlands. This, coupled with the high volume 
addition of water over the spray field was assumed to 
produce a hydraulic mound sufficient to divert the 
background flow around the study area. At lower rates of 
artificial recharge, this would not be the case. 

Aquifer storativity. A rise in the water table 
reflects an increase in the quantity of water stored in the 
aquifer. In a surficial aquifer, this quantity is closely 
approximated by the specific yield of the aquifer. While 
the quantity of stored water ranged widely over the study 
period, its level on June 07, 1988 had returned to its 
January 30, 1988 level, the first day of the study. 


Therefore, this does not enter into this calculation. 
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Final Balance. The calculated values discussed above 
are presented in Table 4-6. This balance indicates an 
unaccounted for loss of 623,000 gallons during a period of 
130 days, or 25 percent of the expected recovered volume. 
However, it should be noted that these quantities are based 


on approximations and average values. 


Table 4-6. Results of System Water Balance Calculations 


Volume pumped from - 1,900,000 gallons 
recovery wells 

Volume injected into 493,000 gallons 
injection wells 

Volume of effective 2,150,000 gallons 
spray on the spray 
field 

Volume of effective 155,000 gallons 
precipitation 

Volume of - 273,000 gallons 


evapotranspiration 


Net volume difference 623,000 gallons 
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Model Selection 


To adequately simulate the Lake Alfred site, a model 
must include the ability to accommodate aquifer 
heterogeneities, irregular flow boundaries, solute 
retardation, and biodegradation to support the bulk of the 
remediation research being conducted on site. Because of 
the well defined contaminant layer across the surface of the 
aquifer, a three-dimensional model would be well suited for 
modeling the site, but their lack of availability and high 
computer memory requirements precluded this possibility. 

BIOPLUME II was selected because it included solute 
retardation, aerobic and anaerobic degradation, and had the 
necessary flexibility to simulate the site. In addition, 
the program was verified, well documented in the literature, 
and supported by the Department of Environmental Science and 
Engineering staff at Rice University. The addition of the 
menu driven preprocessor and its operation on a personal 


computer were two additional attractive features. 


Model Validation 


Before the model could be calibrated, initial 
conditions were required. Once established, the physical 
parameters and contaminant concentrations were entered, and 


the calibration procedure was continued. 
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Initial Conditions 

The BIOPLUME II model, as well as the USGS Solute 
Transport model, require no-flow boundaries to surround the 
modeled region. No-flow boundaries have a transmissivity of 
zero, thus preventing the movement of water or contaminants 
across the cell's boundaries. A constant-head boundary may 
also be assigned to cells, indicating a stable water table 
elevation. Constant-head boundaries, coupled with large 
leakance terms are used to simulate sections of aquifer 
under stable flow conditions. 

The project site was divided into a 200 by 200 foot 
grid system containing 400 square cells, ten feet per side. 
However, 76 of these cells were included in the requisit no- 
flow boundaries. The orientation of the established grid on 
the site map is shown in Figure 4-14. Arbitrary no-flow 
boundaries were established along the northern and southern 
edges of the study area. Since the wash rack proved an 
effective barrier to the northward spread of the plume, a 
boundary was established at its southern edge. Similarly, 
the no-flow boundary along the southern edge was established 
between UF-2M and OHM-2, outside of UF-2M's cone of 
influence. 

A constant head boundary on the western edge of the 
area was established west of RAP-3. The assigned elevation 


was the water level measured in UF-3W. 
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Figure 4-14. Modeled Area of the Lake Alfred Site. 
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On the remaining side, a boundary was established along 
the western edge of the spray field. The boundary location 
was dictated by the line of wells available to generate data 
to set the upgradient constant head boundary. Measured 
water table elevations for RAP-5, RAP-6, RAP-7, and RAP-11 
were entered into their respective cells and values were 
interpolated for the intermediate cells. The additional row 
of cells upgradient are included in the model to minimize 
boundary effects if these four wells are used as injection 
wells. 

By assigning a leakance factor to constant head 
boundaries, a net flux of oxygen and contaminant to or from 
the system may be represented. With a high hydrocarbon 
removal efficiency of the spray aeration system, the 
hydrocarbons not removed by air stripping were assumed to be 
biologically degraded in the upper layer of the soil. 
Therefore, there was no contaminant flux into the system. 
This was a reasonable assumption considering the relatively 
low hydrocarbon contamination levels found in RAP-4, RAP-6, 
and RAP-7. However, the area of the initial gasoline spill 
around RAP-5, had a persistent hydrocarbon concentration. 

To accommodate this contaminant input, the constant head 
boundary in this area was given a leakance coefficient of 
1.0 with contaminant addition, thus allowing the movement of 


contaminant into the area. An influx of oxygen was also 
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used to simulate the movement of oxygenated water into the 
system from the spray field. 

Due to indications of aquifer compression below 
buildings, the transmissivity and thickness of the areas 
underlying these buildings was assumed to be zero. This is 
a reasonable simplification when coupled with the areas of 
relatively high transmissivity surrounding the buildings. 
Paramet Selection 

Killan (1987) determined the effective porosity of the 
aquifer to be 25 percent, and the bulk density to be 1.4 
grams per cubic centimeter. His measured values of aquifer 
transmissivity, presented in Appendix C, were used as the 
original estimates, although they were later adjusted during 
the calibration phase. 

The aquifer thickness map was generated by subtracting 
the surface elevation of the confining layer from the 
average augmented water table elevation in each monitoring 
well, and interpolating between those points over the 
remainder of the area. Thickness values of injection wells 
were disregarded during their period of use. 

The recharge parameter is based on the measured pan 
evaporation and precipitation data obtained from the NWS 
weather station at Lake Alfred. The daily values of these 
measurements for the study period are listed in Appendix G. 
Assuming 10 percent of total rainfall is surface runoff, 


there was an effective rainfall of 13.02 inches during the 
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course of the study. Based on the average annual 
evapotranspiration rate for turfgrass, there was 20.25 
inches of water lost through this process. The net loss of 
7.23 inches of water from grassed areas over the 130 day 
period of the study produces a recharge rate of 5.4 x 0° 
feet per second. Although this value represents a net loss 
from the system, it is given a positive value by the 
numbering convention used in the model. This factor was 
only applied to unimproved surfaces on the site. 

The values of longitudinal dispersivity and the ratio 
of transverse to longitudinal dispersivity were taken from 
Freeberg et al. (1987) as ten feet and 0.1 respectively. 
These values were obtained in a similar calibration study. 

A partition coefficient (Kd) of 1.8 was used for the 
calibration procedure. This value 1s based on laboratory 
experiments conducted by Angley (1987). 

The anaerobic decay coefficient was set at zero. 
Although recent research has found anaerobic degradation of 
alkylbenzenes under restricted conditions, it occurs at a 
slow rate in comparison with aerobic degradation (Wilson et 
al., 1986). Therefore, any anaerobic degradation was 
assumed to be negligible. 

The reaeration decay coefficient was also kept at zero. 
Although there is sure to be some reaeration occurring 
because of the shallow depth of the aquifer, the area of 


highest concentration is effectively sealed under the 








89 
asphalt surface. Infiltration of dissolved oxygen across 
the upgradient constant head boundary is the main source of 
oxygen. 

The process of entering data into the BIOPLUME II 
program is simplified through the use of the menu-driven 
preprocessor. Its use ensures the input file is properly 
formatted, as well as providing the allowable parameter 
ranges. The main loader menu is shown in Figure 4-15. The 
input file may also be constructed or edited with a 
commercial program editor. A sample input file is presented 
in Appendix H, and an edited program output is presented in 
Appendix I. 

Figure 4-16 provides a simplified flow diagram for the 
model. 

Groundwater Flow Calibration 

After adjusting the transmissivity values to minimize 
the summation of errors between the measured water table 
elevations and the predicted water table elevations in the 
eight monitoring wells under the three flow conditions, the 
minimum error obtained was 6.108 feet. A breakdown of this 
value is provided in Table 4-7. 

Sources of error include geologic heterogeneities in 
the aquifer and variations in the hydraulic gradient not 
accounted for in the model. The model also uses nodes 
centrally located in the finite difference cells. Due to 


the grid-size limitation of the model, there is a maximum of 
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Loader 


Main Menu 


Oe Edit file name 

De Edit card 1 (Title) 

oe Edit card 2 (Grid/timing parameters ) 

a Edit card 3 (Grid/timing parameters ) 

Ds Edit card 4 (Reaction parameters) 

6. Edit data set 1 (Observation wells) 

72 Edit data set 2 (Pump/Inject wells) 

8. Edit data set 3 (Transmissivity map) 

On, Edit data set 4 (Thickness map) 

nO, Edit data set 5 (Recharge map) 
Agdie Edit data set 6 (Nodeid map) 

te Edit data set 7 (Nodeid code definitions) 
ee. Edit data set 8 (Water table elevations) 
14. Edit data set 9 (Initial hydrocarbon conc. ) 
allen Edit data set 10 (Initial oxygen conc.) 
16. Edit data set 11 (Pumping periods) 

17. Write data to file 
1385 Quit 


Enter the number of your choice (1..18) 


Figure 4-15. BIOPLUME II Preprocessor's Main Menu 





GENERATE UNIFORM 
DISTRIBUTION OF TRACER 


DETERMINE LENGTH OF TINE 
INCREMENT FOR EXPLICIT 
CALCULATIONS 


MOVE CONTAMIWAIT 
PARTICLES 


GEBERATE OR RENOVE 
PARTICLES AT 
| APPROPRIATE BOUNDARIES 


COMPUTE COMWTANINANT 


| COMCENTRATION IW CELLS 
AND AT NODES 


COMPUTE MASS BALANCE 


SUPERFOSE CONTAMIFANT 
ABD OXYGEN PARTICLES 


ADJUST MASS BALANCE 


ERD OF TIME 
STEP? 


Figure 4-16. 
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Simplified Flowchart of BIOPLUME ITI 
(Source: Rifai et al., 


1987, p. 2-9). 
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Table 4-7. Minimized Groundwater Calibration Error 


Condition 


Wel Static RAP-5 and RAP-7 Peon aaa i. 
RAP-4 02265 -0.067 67037 
RAP-12 -0.218 -0.649 —On2 91 
Mai Onis 0.083 Oma o 
RAP-13 O2577 O72o7 0.139 
UF-—15 0.252 -0.024 0.184 
P-5 =O2 15 -0.451 -0.398 
OHM- 3 -0.070 -0.036 02633 
OHM-4 Onl33 -0.033 0.534 


Minimized Total Error 


+ [P, - M| 


6.108 feet 
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seven feet between the actual location of the monitoring 
well in the field and the mathematical representation of the 
well in the model. Considering these factors, the 
variability of the water table elevation, and the similar 
results obtained under three diverse flow conditions, the 
model is considered adequately calibrated. 

Solute Transport Calibration. 

Meta-xylene and para-xylene are present in sufficient 
quantities in the contaminated groundwater to serve as 
indicators of the plume movement; therefore, they were used 
to calibrate the solute transport portion of the model. The 
sum of both isomers was used since the gas chromatography 
method used for analysis failed to differentiate between 
them. Although the same argument could be made for ortho- 
xylene, the meta- and para- isomers were selected because of 
their faster degradation rate (Angley, 1987). This tracer 
method will be effective under conditions without 
degradation, but where degradation occurs, the oxygen will 
be consumed in the degradation process of all organics, not 
just the tracer organics. 

The unknown contaminant distribution or concentration 
beneath the asphalt roadway, the limited database available, 
and the instability of the database that was available gave 
additional sources of error. 

BIOPLUME II does not effectively model this site for 


solute transport. The extreme vertical variations of 
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contaminant concentration involved, make it unrealistic to 
predict the contaminant concentration in the groundwater. 
The highly contaminated regions serve as an intermittent 
source of contaminant whenever the saturated zone enters 
these areas. Preliminary calculations of the model, and the 
elevated hydrocarbon concentrations observed after the high 
volume injections in P-6 and P-7 indicate combined 
concentrations of 150,000 to 250,000 parts per billion of 
meta-xylene and para-xylene in the area beneath the asphalt 


surface. 





CHAPTER V 
SUMMARY AND CONCLUSIONS 

The objectives of this study have been met. A semi- 
quantitative evaluation of the vertical distribution of 
sorbed hydrocarbons in the soil was conducted, confirming 
that high concentrations of solute are localized in areas 
previously contaminated with free product. The vertical 
stability of the sorbed contaminants is further supported by 
the large contaminant concentration swings in the 
groundwater which correlate with the fluctuations in the 
water table elevation. 

A mass balance of the water flowing through the system 
was performed, using pan evaporation and precipitation 
values, pumped water volumes, and sprinkler application 
efficiencies. There was an unaccounted for loss of 25 
percent of the estimated recoverable water, a reasonable 
margin of error considering the estimations involved in the 
process. This justified considering the study area as a 
closed system for the purposes of modeling. The most likely 
avenue of unaccounted water loss is from the north face of 
the spray field toward building 14. If the assumptions of 


total hydrocarbon removal through air stripping and 
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biological degradation in the surface layer of soil are 
valid, there will be no adverse impact. 

BIOPLUME II was selected as the modeling program based 
on its adaptability to variable site conditions, thorough 
documentation of the program and its predecessor, the USGS 
Solute Transport Model, the program's availability and 
institutional support, and the user friendly software. 

The program was calibrated to the site conditions using 
a trial and error method to minimize the difference between 
predicted and actual values. Despite the physical 
obstructions and heterogeneities within the aquifer, the 
program was able to readily simulate the groundwater flow 
conditions; however, the solute transport model was not able 
to be adequately calibrated. To account for the highly 
contaminated region at the surface of the elevated water 
table, intermittent contaminant sources would be required, 
but the concentration at the source would also depend on the 
water table elevation and would vary over time. While this 
aids the overall model, it reduces the model's ability to 
realistically predict the dissolved contaminant movement. 
This inability to predict reasonable transport will severely 
restrict its usefulness. To obtain a more sensitive model, 
a three-dimensional model would be required. 

The database available for model calibration at 
elevated water table conditions was very limited. 


Fluctuations in the water table significantly affect the 
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contaminant concentrations. The transmissivity values of 
the area remain constant. 

The largest area of contamination remains under the 
road surface between buildings 12 and 31. Flushing the 
aquifer via fresh water injection proved effective at 
removing large quantities of contaminants. Because of the 
low cost and high effectiveness, flushing under the road 
surface, with recovery wells operating, should be an 
effective remedy on high contaminant concentrations. Follow 
on work with bioremediation should be much more effective 
with lower concentrations remaining after flushing. 

Development of a more extensive database could improve 
the calibration of the model. In particular, improved 
information on the biological parameters and quantification 
of the distribution and concentration of contaminants under 


the road surface would improve the reliability of the model. 








APPENDIX A 


RECORD OF WATER TABLE ELEVATIONS 
IN MONITORING WELLS 





Table A-1. Water Table Elevations 


Elevation (feet, a.m.s.1.) on Specified Date 


-_— 30 Jan 88 O2 Feb 88 06 Feb 88 O09 Feb 88 13 Feb 88 
M2 139.28 136.40 136.45 136.30 135295 
M-2 i384 .72 135852 Lsi5e65 13552 W3i5e2 7 
M=-3 134.81 £35.63 135°. 73 135-63 135.39 
OHM-1 13358 T3238 133266 1332.90 Ali e) stexe) 
OHM-2 igsy 2s 23523 ats 33). 05 i3i2. 90 
OHM- 3 Anes Cis fal 134.04 134519 134.14 134.04 
OHM-4 Leeeo9 iesroZ 133299 133.44 133.34 
P=5 132 #82 13). 62 135.74 Tei) eye) 135.34 
2-6 iin 08 136500 13620 13 5n2 5 153'5)..6 3 
Be 7 iso. 19 16. 29 136m39 iGo.26 135.66 
RAP-—2 133% aL L3sm06 see 1 133.74 13269 
RAP-4 136.06 U3? 28 ghey 5 ore Peie23 136. 76 
RAP-5 136.70 1432.03 144.13 144.13 144.13 
RAP-6 ian 07 iso. 34 sens 2 338.29 37.77 
RAP-7 136.34 140.44 142.14 142.14 142.14 
RAP-8 =o 145.52 145.52 144.12 fa2.37 
RAP=9 dg. 52 35.97 P38..09 Se. a7 138.35 
RAP-10 PES) 7 4a igews1 138.66 1337 6 e38-11 
RAP-11 136 +80 Leena 7 138332 138.42 Seis We) 
RAP-12 TS6e02 Leis 7) LSa365 Let Pons? 
RAP-13 134.8 T36)..05 ie4. 13 136.00 USS 55 
UP-1E —— = == UZG202 gs. 72 
UF-3W L325 9 132.64 132.74 132.47 i32.09 


Ne 
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Table A-i--continued. 


Elevation (feet, a.m.s.1.) on Specified D 


Os 6p Pep 83 20 Feb 8&6 23 Feb 88 27 Feb 88 01 Mar 88 
m= 1 136.08 a ey) 5 1518) Nese SIS) 135.43 Poon 2o 
M-2 135720 Woo 22 134.98 134.75 134.62 
m3 35.33 Ps'5.36 iS 5:12 134.86 34.73 
OHM-1 dv3i3. 61S 133.48 P3S.cs 133.30 HEGRE (Uke) 
OHM-2 12.33 133.00 ae W327 5 132.60 
OHM-3 3 oa.99 134.44 133.94 133.54 133.64 
OHM - 4 133.34 pC ee ayy 133.24 133.04 PS 2eo4 
Se HS 5.37 i529 135.06 134.84 134.72 
Eo 135.68 $35).53 1361 3 59.1.3 134.93 
eal tea. 91 135.74 135.49 12s © 139.06 
RAP-2 133.61 33.09 133.83 133.24 133224 
RAP-4 136.96 136.61 136.36 136.28 136.08 
RAP-5 144.13 144.13 144.13 144.13 142.43 
RAP-6 33.14 phy eis L320 136.52 136212 
RAP—7 141.36 im 14 142.14 142.14 141.44 
Par -8 122.52 141.42 ao 141.62 al 27 
RAP-3 136.65 USS 3s7 137.67 138.07 137 215 
RAP-10 138. 54 35:09 138.41 L37.o1 137.44 
RAP -11 MSS 0 37310 P36. 79 37.40 137.02 
RAP-12 JS) TF gal eG. 92 Psi6562 1S Spans Te) Sil LZ 
RAP-13 135.60 135.43 135).10 P3o..05 134.73 
UF-1E 35.07 PS5562 £35.35 P35.22 134.87 


UF-3W 132.409 132.34 132.34 132.22 Dea Geet 
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Table A-1--continued. 


Elevation (feet, a.m.s.1.) on Specified Date 


. Q2 Mar 88 05 Mar 88 08 Mar 88 12 Mar 88 15 Mar 88 
M-1 3.5 45 135.40 135.61 135276 136228 
M-2 135407 134.70 134.87 135.00 136.30 
M-3 ioe L3 134.81 136.00 ACs i! 134.58 
OHM-1 ig 3.35 NBS... 23 132.102 33.75 134.28 
OHM-2 32 260 132055 132.80 132278 133.18 
OHM-3 3 264 133.54 a3 09 133.76 133796 
OHM-4 313 40:7 132.94 133.08 133.14 133739 
P-5 225 oT 134.77 134.97 135.08 1357.33 
P-6 S559 135.03 35. 23 Alc usie 61S} 13567 
P-7 135.74 iso 19 135.68 35.49 135.83 
RAP-2 23526 33526 133641 133.46 133.68 
RAP-4 PG. 93 SOc 3 136.46 6. OL 136.81 
RAP-5 142.80 m2 .55 -- 142.76 143.06 
RAP-6 we 07 LS'7 S34 137.64 DS PeTS 37792 
RAP-7 141.79 141.89 = 142.04 142.14 
RAP-8 aa. 72 141.62 = 141.82 141.97 
RAP-9 138 282 137 287 138.07 38.27 HUGE eee 
RAP-10 iS 8516 Ley 261 779 137799 138.24 
RAP-11 PSS a2 137 222 137.39 Sto 136.75 
RAP-12 M37 012 136 532 136.47 136662 136.90 
RAP-13 135.50 134 235 rS'505 Sel 6 BES ee Li) 
UF-1E 1350537 134.97 PS 5q 22 eS oee 7 CBS fue 


UF-3W i 131299 132-16 13299 32501 
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Table A-1--continued. 


Elevation (feet, a.m.s.].) on Specified Date 


a 19 Mar 88 22 Mar 88 26 Mar 88 30 Mar 88 0Q2 Apr 88 
Me 1 136. 58 SE 1357.92 135.80 135.78 
M-2 15e72 -- 135205 135.27 lg Se) 
m3 P35". 82 a= 135.18 135.33 134.43 
OHM-1 133.60 134.76 13145465 134.78 134.64 
OHM- 2 res 740 3320 132793 133.45 133719 
OHM-3 134.24 134.14 33295 134.12 133.84 
OHM- 4 33. 59 133.49 ss. 22 vse. 79 Lee 4 2B 
P-5 13S", 82 135.46 135°. 12 35.32 135.16 
Pat 136.14 13565 135.43 135258 135.45 
re] 136. 39 136. 79 135.58 1357-69 135.67 
RAP-2 33. 96 13395 P3358 esi} 133.256 
RAP-4 137.46 137.48 136.66 136.43 136259 
RAP-5 143.41 139.30 139..28 141.13 142.75 
RAP-6 138.69 een To 137 sl 37521 137.62 
RAP-7 142.14 141.64 141.66 137.74 1A 58 
RAP-8 42.25 141.90 142.32 142.42 142.37 
RAP-9 139% 27 13. 52 138.41 136.52 138.27 
RAP-10 #38. 916 138. 21 138.10 137331 137799 
RAP-11 138. 52 137.82 137.67 P7707 137.09 
RAP-12 137.59 136.87 136.72 136.52 136.72 
RAP-13 136.10 135.48 130.26 35.0 130.00 
UF-1E 136.31 es 567 153:5513.7 1S S57 135259 


UF-3W inst. 42 132.49 32.20 32249 32.49 
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Table A-1--continued. 


Elevation (feet, a.m.s.l.) on Specified Date 


a O05 Apr 88 09 Apr 88 17 Apr 88 24 Apr 88 29 Apr 88 
M-1 136.04 139. 08 135.43 #30703 Poon a7 
mM—2 135.34 134.90 134.79 134.55 is5269 
M-3 135.43 30.02 134.92 134.68 135.64 
OHM-1 134.70 134250 134.00 133295 34.27 
OHM-2 3.10 S203 132.64 SLs IE, 32.62 
OHM-3 134.20 23576 is267 33.07 134.00 
OHM-4 isae'S3 133.16 SSCs 33.38 Was.o1 
e='5 135.45 is 29'S 134.86 134.65 135.94 
Eo nS. 78 135224 135.08 134261 PsSeco 
od sd 99 iS L TP3o.22 134.89 138590 
RAP-2 133.86 1373.49 P33e39 se o2 133261 
RAP-4 136.86 ioOn so PSG 21 TS ieas 4/16) Voy soe 
RAP-5 ia . 93 4a2el2 144.13 136.14 37390 
RAP-6 Leh Galil cls Wyre 24 US} 7/ ail £36.51 1372782 
RAP-7 141.56 141.42 142.14 Hele ayek 137.00 
RAP-8 142.18 141.09 141.85 142.40 141.57 
Rae =9 138.67 dea 7 97 1S. 76 136.69 138.37 
RAP-10 ions 7 iS. 66 137.49 136.48 138.14 
RRP-11 ey 294 PST e2 6 137.12 136.20 137.84 
RAP-12 S712 136.45 sG2 2 i310. 67 ots 39 
RAP-13 do. 66 elias 134.88 134.49 137.00 
UF-1E 35.85 £35224 134.97 134.72 136.16 


UF-3W woz. 31 eZ. 2 131.74 PS Ps07 30,907 
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Table A-1--continued. 


Elevation (feet, a.m.s.1.) on Specified Date 


Bes O02 May 88 O06 May 88 14 May 88 20 May 88 0QO7 Jun 88 
M-1 1364.08 137.07 S572 35.07 134.83 
MZ 1397.98 136.35 135.97 134.94 134.42 
m-3 135.94 136.40 s4. 59 135). 10 134.58 
OHM-1 134.67 134.56 134.52 134.29 134..30 
OHM-2 s2.9] NG PAB Js) eae 17 132.48 132.44 
OHM-3 134s 31 134.43 134.05 MEN nS) 13 31 O10 
OHM-4 i3s85 133266 3539 S309 132.87 
E> so. 21 136.40 13a. 30 134.99 134.49 
r=6 desi o 7 P3920 39516 35.15 134.63 
el 4021 140.14 135.74 135.23 134.74 
RAP-2 33289 133.99 133g.383 134.01 133047 
RAP-4 39753 137.66 136.46 136.30 PSS 201 
RAP-5 ise 32 1367.11 137.05 136.89 13 Gn 
RAP-6 138.14 t8.29 137.14 Ne 6 Si) USGn22 
RAP-7 7 24 V37.37 136.65 36.7) seis 5 12/7 
RAP-8 == 139245 138.20 138262 137.47 
RAP-9 138.80 138.64 137.50 137.594 136.63 
RAP-10 38. 59 133.41 27 eae) 137329 136.40 
RAP —-11 S023 138223 137 202 SO od SG. L2 
RAP-12 o7 94 7.78 HbSiO- 92 136722 ee 
RAP-13 2S 237.05 35.0 134.89 134.35 
OE-1E ES66.2 136), a1 135.40 134.94 134.62 


UF-3W ESie8 2 132.103 oZ202 131.16 Dole72 








APPENDIX B 


PUMPING RECORD OF THE GROUNDWATER RECYCLING SYSTEM 





Bable B=-1. 


Date 


30 
31 
01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
lial 
a2 
a3 
14 
15 
16 
a7 
18 
19 
20 


2a 


Jan 
Jan 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 
Feb 


Feb 


Project 
(1988) Day 


W271 
728 
29 
730 
731 
732 
ues 
734 
25 
736 
137 
738 
39 
740 
741 


742 


744 
745 
746 
747 
748 


749 


Pumping Record of the Groundwater Recycling 


Normalized Daily Pumped Volume (gallons) 


RAP-1 


6631 
7015 
8473 
S136 
7200 
7200 
£3931 
10016 
$972 
8523 
8077 
6500 
5700 
6656 
6969 
707 3 
TOTS 
6678 
3269 
2718 
0 

0 

0 


UF-2M 


5778 
7298 
1298 
TN eas 
8000 
8000 
8776 
8789 
8793 
8793 
8358 
8377 
8413 
8464 
8772 
8863 


106 


RAP-3 


1051 


Spray 
18637 
17631 
11246 
0951 
12000 
12000 
9683 
2459 
Bog 
13200 
13200 
13000 
13000 
13000 
14213 
14493 
14493 
11416 
POs 
10398 
Ls 10 
P67 2 


11779 








Wiyitersies ecosts oes 


if 
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Normalized Daily Pumped Volume (gallons) 


@epie B-1 -- continued. 
Date Project 
(1988) Day RAP-1 
22 Feb 750 1142 
23 Feb (opt 8687 
24 Feb 752 12825 
25 Feb 158 12825 
26 Feb 754 12825 
27 Feb TOS 4183 
28 Feb 756 1140 
29 Feb mT 1140 
01 Mar Tas 7116 
02 Mar 759 8037 
03 Mar 760 6506 
04 Mar 761 6733 
05 Mar 762 7323 
06 Mar 763 7416 
07 Mar 764 8119 
08 Mar 765 8531 
09 Mar 766 8623 
10 Mar 767 9243 
11 Mar 768 9371 
12 Mar 769 7972 
13 Mar 770 7649 
14 Mar 1 7959 
15 Mar 772 7963 


UF =2M 


8770 
8624 
8533 
8533 
8533 
8735 
7987 
7987 
1233 
7364 
7350 
7297 
7806 
7886 
7920 
7719 
7623 
8020 
8019 
7417 
208 
7617 


6774 


RAPS 3. 


490 
346 
582 
582 
582 
785 
Si5i9 
855 
451 
378 
418 
429 
1102 
1207 
1658 
2049 
2167 
2205 
1325 
1135 
L092 
1068 
1345 


spray 
Tr399 
11427 
11446 
11446 
11446 
11341 
11305 
Piss 
11242 
11341 
12i22 
11051 
nes] 
10989 
11830 
14129 
15060 
£3236 
13161 
13100 
13085 
2379 


2n53 








Table B-1 -- continued. 
Date Project 
eigss) Day  _RAP=1_ 
16 Mar 173 7925 
17 Mar 774 8246 
18 Mar 175 8784 
19 Mar 776 6565 
20 Mar 1 6186 
21 Mar 778 6186 
22 Mar 779 6106 
23 Mar 780 6045 
24 Mar 781 8170 
25 Mar 782 10581 
26 Mar 783 7026 
27 Mar 784 5969 
28 Mar 785 5764 
29 Mar 786 6023 
30 Mar 787 6120 
31 Mar 788 6490 
01 Apr 789 8953 
02 Apr 790 5757 
03 Apr 791 4968 
04 Apr 792 5122 
05 Apr 793 5965 
06 Apr 794 6900 
07 Apr 795 8100 
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Normalized Daily Pumped Volume (gallons) 


UE=2M 


5545 
5032 
4828 
4838 
4840 
4840 
4670 
4192 
3968 
3710 
2370 
1979 

200 


0 


5200 
4600 


RAP-3 
2837 
3096 
Soe 
1208 

836 
836 
816 
680 
1643 


1772 


807 


882 


679 
727 
887 
1720 
1670 


Spray 
14000 
129377 
12864 
12792 
12780 
12780 
12782 
12700 
12688 
12642 
12844 
12912 
12410 

9255 

6576 
11439 
10814 
10551 
10477 
10452 
10703 
10280 
10690 





Table B-1 -- continued. 
Date Project 
(1988) Day RAP-1 
08 Apr 796 8621 
09 Apr 7197 Sino2 
10 Apr 798 4542 
11 Apr 799 4500 
12 Apr 800 5300 
13 Apr 801 5200 
14 Apr 802 TI2 
u5 Apr 803 7932 
16 Apr 804 7932 
17 Apr 805 5415 
18 Apr 806 4700 
19 Apr 807 6400 
20 Apr 807 8400 
21 Apr 809 8300 
22 Apr 810 8183 
23 Apr 811 8183 
24 Apr 812 6135 
25 Apr 813 4500 
26 Apr 814 5400 
27 Apr 815 7100 
28 Apr 816 9147 
29 Apr 817 8783 
30 Apr 818 8685 


4891 
4518 
4394 
4400 
4600 
4500 
4940 
4940 
4940 
4582 
5100 
3700 

0 


4400 
6000 
7200 
8000 
fie 
7638 


1477 
422 
v1 
1500 
1320 
820 
19ar3 
iors 
LGars 
1044 
690 
1770 
3190 
3160 
2947 
2947 
2026 
2740 
3300 
3680 
3724 
2655 
2307 


i} 


Normalized Daily Pumped Volume (gallons) 
UF-2M 


Spray 
10122 
iesz0 
10398 
10380 
9960 
10330 
10144 
10144 
10144 
10107 
9920 
10120 
10640 
9310 
8919 
8919 
8711 
8550 
9340 
9270 
9280 
9425 


9468 





Table B-1 -- continued. 
Date Project 

(1988) Day RAP-1 
01 May 819 8685 
02 May 820 9245 
03 May S21 9203 
04 May 822 9752 
05 May 823 9787 
06 May 824 8698 
07 May 825 8698 
08 May 826 8698 
09 May 827 8854 
10 May 827 9530 
11 May 829 9835 
12 May 830 9419 
13 May 831 9030 
14 May 832 8201 
15 May 833 7983 
16 May 834 8068 
17 May 835 7905 
18 May 836 8034 
19 May 837 8469 
20 May 838 7650 
21 May 839 7650 
22 May 840 7650 
23 May 841 7595 


7638 
8029 
7851 
7928 
7958 
7347 
T3847 
7347 
6930 
6982 
6998 
6822 
6789 
6470 
6380 
6343 
6152 
6118 
6099 
5948 
5948 
5948 


5874 


RAP-3 


2337 
4622 
3360 
3646 
3867 
2641 
2641 
2641 
3068 
3354 
3700 
2608 
2736 
2397 
2305 
2136 
2114 
1889 
1887 
2025 
2025 
2023 


Ale taes) 
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Normalized Daily Pumped Volume (gallons) 
UF-2M 


spray 
9468 
9232 
9331 
9293 
9122 
9215 
9215 
9215 
9156 
4189 
1340 
0 

0 
5184 
9570 
9478 
9410 
8787 
9679 
9491 
9491 
9491 
9479 
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Pable B-1 -- continued. 

Date Project Normalized Daily Pumped Volume (gallons) 
(1988) Day RAP-1 UF-2M RAP-3 spray 
24 May 842 8494 6245 2736 9439 
25 May 843 9548 6722 4207 909 
26 May 844 9325 6386 3828 0 
27 May 845 9360 6172 3360 0 
28 May 846 9360 6172 3360 0 
29 May 847 9360 6172 3360 0 
30 May 848 9360 GliZ 3360 0 
31 May 849 8023 5965 ee) 0 
01 Jun 850 8646 6010 1610 0 
02 Jun 851 8646 5979 1620 0 
03 Jun 852 8646 5866 1978 0 
04 Jun 853 8646 5866 1978 0 
05 Jun 854 8646 5866 1978 0 
06 Jun 855 8646 5oe5 2190 0 
07 Jun 856 7264 5500 1460 6672 


Table B-2. 

Date Project 
(1988) Day 
30 Jan T27 
31 Jan 728 
01 Feb 729 
02 Feb T1350 
03 Feb 731 
04 Feb nS 2 
05 Feb 733 
06 Feb 734 
07 Feb 135 
08 Feb 736 
09 Feb 737 
10 Feb 738 
11 Feb 739 
12 Feb 740 
13 Feb 741 
14 Feb 742 
15 Feb 743 
16 Feb 744 
17 Feb 745 
18 Feb 746 
19 Feb 747 
20 Feb 748 
21 Feb 749 


Injection Record in the Well System 


pil 


Normalized Daily Injected Volume (gallons) 


REP- 5 


ay. 
6863 
6477 
6578 
7303 
7303 
7441 
3836 
3004 
5894 
3430 
6240 
6640 
5185 
3067 
2510 
2010 
2454 
2464 
2470 
2677 
1893 


1699 


RAP-7 


jal Tia 
6863 
6477 
6578 
7303 
7303 
7441 
3836 
3004 
5894 
3430 
6240 
6640 
5iS5 
3067 
2510 
2510 
2454 
2464 
2470 
Zio 7 
S23 
1622 


RAP-8 
8417 
12243 
11085 
6663 
4800 
4800 
4909 
2514 
1962 
1708 
985 
800 
700 
610 
714 
741 
741 
697 
651 
647 
726 
736 
738 


eae @f fo of © 


OOO On © an anc 


en 


p- 


0206CODWUmUmCOCOWUCOTCUCOClCOUlCOUlCOCClCOTT CTC Cc 5S— 9 5(—€@«0«' Fs BD 


(om) 


oO olUlmUOmUmlmlCUCODUCO 


j 
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Table B-2 -- continued. 
Date Project 
(1988) Day RP= 5 
22 Feb 750 1726 
23 Feb 751 1825 
24 Feb 752 1877 
25 Feb 753 1877 
26 Feb 754 1877 
27 Feb 70 982 
28 Feb 756 707 
29 Feb 757 707 
01 Mar 758 703 
02 Mar 759 659 
02 Mar 760 659 
04 Mar 761 659 
05 Mar 762 616 
06 Mar 763 609 
07 Mar 764 609 
08 Mar 765 602 
09 Mar 766 598 
10 Mar 767 598 
11 Mar 768 598 
12 Mar 769 597 
i Mar 770 a7 
14 Mar 771 O97 
15 Mar 772 618 


1726 
1825 
Leo? 
1877 
1877 
982 
707 
707 
703 
659 
659 
659 
616 
609 
609 


602 


SS) 
598 
597 
oa7 
597 


618 


RAP-8 


Use 
682 
662 
662 
662 
630 
621 
621 
611 
525 
525 
525 
482 
475 
475 
444 
430 
430 
430 
409 
404 
404 


413 


P-6 


0 


0 


oO Oo Oo 


(oo) 
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Normalized Daily Injected Volume (gallons) 
RAP-7 


Pe? 
0 


0 


(=) 


oO Oo oOo 


seer eg eg 








‘ ;crogF 
; RR ' | Bak 
, ’ ' 


Table B-2 -- continued. 
Date Project 
(1988) Day RAP-5 
16 Mar 773 632 
17 Mar 774 632 
18 Mar 775 Gaz 
19 Mar 776 464 
20 Mar Ciel 423 
21 Mar 778 423 
22 Mar 779 407 
Za Mar 780 so7 
24 Mar 781 397 
25 Mar 782 397 
26 Mar 783 592 
27 Mar 784 644 
28 Mar 785 644 
29 Mar 786 644 
30 Mar 787 641 
31 Mar 788 640 
01 Apr 789 640 
O02 Apr 790 660 
03 Apr 791 602 
04 Apr 792 602 
05 Apr 793 575 
06 Apr 794 562 
07 Apr 795 562 


RAP-7 
632 
632 
632 
464 
423 
423 
407 
oo7 
S27 


307 


644 
644 
644 
641 
640 
640 
660 
602 


602 


RAP-38 
418 
418 
418 
291 
259 
209 
409 
499 
499 
499 
e)S))0) 
354 
354 
354 
344 
338 
338 
347 
349 
349 
334 
326 


326 


Peo 


0 


0 
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Normalized Daily Injected Volume (gallons) 


P-7 
0 


0 





Table B-2 -- continued. 
Date Project 

(1988) __Day RAP-5 
08 Apr 796 562 
09 Apr To 1208 
10 Apr 798 1461 
11 Apr 799 1461 
12 Apr 800 1461 
13 Apr 801 1259 
14 Apr 802 1236 
iS) Apr 803 1236 
16 Apr 804 1236 
17 Apr 805 571 
18 Apr 806 447 
19 Apr 807 447 
20 Apr 808 447 
21 Apr 809 447 
22 Apr 810 447 
23 Apr 811 447 
24 Apr 812 0 
25 Apr aus 0 
26 Apr 814 0 
27 Apr 815 0 
28 Apr 816 0 
29 Apr 817 0 
30 Apr 818 0 


262 
1208 
1461 
1461 
1461 
1259 
1236 
1236 
1236 

S74 

447 

447 

447 

447 

447 


447 


RAP-8 
326 
319 
316 
316 
316 
316 
316 
316 
316 
272 
264 
264 
264 
264 
264 
264 
264 
264 
264 
264 


264 


PES 


0 


6852 
6852 
6852 
7017 


7066 


shall) 


Normalized Daily Injected Volume (gallons) 
RAP-7 


Bead) 
0 
0 


7449 
7449 
7449 
7576 


7614 
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Table B-2 -- continued. 


Date 


01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
aT 
we 


13 


May 
May 
May 
May 
May 
May 
May 
May 
May 
May 
May 
May 


May 


Project 
(1988) Day 


819 
820 
821 
822 
823 
824 
825 
826 
827 
828 
829 
830 


831 


0 


0 


RAP-7 
0 
0 


RAP-8 
0 
0 


oreo 
7066 
6203 
5912 
3929 
Deo2 
4476 
4476 
4476 
4281 
4189 
4271 
4269 


4242 
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Normalized Daily Injected Volume (gallons) 
RAP-5 


a del 
7614 
6738 
6362 
63390 
6305 
0 


0 





APPENDIX C 


HYDRAULIC CHARACTERIZATION OF CONTAMINATION SITE 
LAKE ALFRED, FLORIDA 





Table C-1. Summary of Estimates of Hydraulic Parameters for 
the IFAS-CREC Project Aquifer. 


Relevant Hydraulic Specific Seepage 
Area of Conductivity Transmissivity Yield Velocity 
Application (fpd) (qpd/ft) (%) (fpd) 
Within 40 ft 
of UF-2M 160 4,700 B0) 18 
Within 80 ft 
of RAP-1 310 12,400 110) 10 
Between RAP-9 
and RAP-10 30 400 A 16 
Between RAP-10 
and RAP-11 26 580 21 10 
Between RAP-11 
sed P27 Sik 580 ale 6 
Between P-7 
and UF-2M 60 810 30 10 
Near P-5 11 220 Bal Zz 
Near P-6 7 130 21 1 
Near P-7 : 8 140 21 2 
Near UF-3W 4 170 30 <1 


Source: Killan, 1987. 
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APPENDIX D 


GAS CHROMATOGRAPHY CALIBRATION CURVES 
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APPENDIX E 


HEXANE EXTRACTION RESULTS OF SOIL CORE SAMPLES 
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Table E-5. Extraction Results from Soil Core 5, North of 


UF-2M. 
Measured Concentration 
Analyte (nanograms hydrocarbon/gram_ soil) 
Benzene BDL (550) BDL (378) BDL (370) 
Toluene BDL (431) BDL (296) BDL (290) 
Ethylbenzene 1170 BDL (377) BDL (369) 
Meta-,Para- 1130 BDL (326) BDL (319) 
xylene 
Ortho-xylene 1230 BDL (292) ND (0) 
Isopropyl- 1080 BDL (295) BDL (289) 
benzene 
Propylbenzene 1430 BDL (331) BDL (325) 
3,4-Ethyltoluene 1110 BDL (321) BDL (315) 
1,3,5-Trimethyl- BDL (388) 932 BDE (262) 
benzene 
2-Ethyltoluene BDL (401) BDL (276) ND (0) 
1,2,4-Trimethyl- BDL (433) BDL (298) BDL (292) 
benzene 
1,2,3-Trimethyl- ND (0) ND (0) ND (0) 
benzene 
Total 9340 4120 2830 


Mean Elevation of 
Sample (feet above 13407 134.0 3364 
mean sea level) 





Table E-5 -- continued. 
Measured Concentration 
Analyte (nanograms hydrocarbon/gram soil) 
Benzene BDL (488) 62900 152 
Toluene 2790 338000 3120 
Ethylbenzene 5720 186000 5180 
Meta-,Para- 18800 571000 15000 
xylene 
Ortho-xylene 8180 234000 7000 
Isopropyl- 1040 22600 1080 
benzene 
Propylbenzene 3790 87300 2570 
3,4-Ethyltoluene 15400 380000 9380 
1,3,5-Trimethyl- 5980 143000 3520 
benzene 
2-Ethyltoluene 7470 139000 3910 
1,2,4-Trimethyl1- 18400 420000 10800 
benzene 
1,2,3-Trimethyl- 4780 93900 2720 
benzene 
Total 92800 2680000 65000 


Mean Elevation of 
Sample (feet above US Peet Uist gts} 12220 
mean sea level) 
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Table E-5 -- continued. 
Measured Concentration 
Analyte (nanograms hydrocarbon/gram soil) 
Benzene BDL (226) 
Toluene BDL (177) 
Ethylbenzene 13S 
Meta-, Para- ASSO 
xylene 
Ortho-xylene 729 
Isopropyl- BDL (176) 
benzene 
Propylbenzene BDL (198) 
3,4-Ethyltoluene 878 
1,3,5-Trimethyl- BDL (159) 
benzene 
2-Ethyltoluene 428 
1,2,4-Trimethyl- 1260 
benzene 
1,2,3-Trimethyl- 5h 7/ 
benzene 
Total 7050 


Mean Elevation of 
Sample (feet above 130.4 
mean sea level) 





APPENDIX F 


RESULTS OF EXTRACTION RECOVERY STUDIES 
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APPENDIX G 


CLIMATOLOGICAL DATA FOR LAKE ALFRED, FLORIDA 
30 JANUARY TO 07 JUNE 1988 





Table G-1. Daily Measurements of Pan Evaporation at Lake 
Alfred, Florida, 30 January to 07 June 1988. 


Day of Evaporation (inches } 
Month January February March April May June 





a Orie OmZz0 0.24 0.04 Ona 
2 O..13 0.26 0.28 0.44 Or au) 
3 0.14 Gri? Oe C35 0.41 
4 OFS ale alae 0.16 Cieze OS) 
2) 0.14 (6) 5 aipk Ors 0.41 0.24 
6 0.07 0.04 0.37 0.43 0.14 
i 0.16 0.20 0.41 Orsl OL 1S 
8 0.05 0.26 0.30 07316 

9 0.20 Omat Crag O30 

10 0.36 Os2 O25 

otal Oeil Cre? O228 Orc 

EZ Od 0.38 Ont Ones) 

3 Geis Ol ale) Oni O.12 

14 Grit 0.43 0.25 Om2y, 

5 0.10 0.28 0.34 OrZi 

16 0.14 O22 0.24 Ome 

i) Ons Slats) O220 0.24 

18 Out? 0.30 O33 

BO 0.16 Oie2/ O25 0.31 

20 0.30 OR22 0.25 0.40 

2a. Orit Orn? 0733 0.43 

Ze 0.27 0.25 OS Pa, 0.28 

23 Oma 0.24 0229 

24 G.27 Os 583 0.30 0.54 

25 Oe22 0726 Once Oni2 

26 Oral 0.30 0.28 

aH 0.24 0.19 0.39 Oz2Z6 

28 O20 022s 0.42 0.35 

20 0.19 0.26 Oe2Z3 0239 

30 Oi2 ero Be al Ome7 

oa G16 0.30 0.34 
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Table G-2. 


Day of 


Month 


Daily Measurements of Precipitation at Lake 
Alfred, Florida, 


30 January to 07 June 1988. 
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Precipitation (inches) 
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APPENDIX H 


SAMPLE INPUT FILE FOR BIOPLUME II 





run-i9 Based on 06 May Data, with 30 Apr - 06 May pumping data 


1 1 20 208000 1 ia 0 100 5 9 3 0 0 0 0 tt) 1 
02 .0010 .250 2.0 .00 .00 O. 10.0 10.0 .1001.000 1.00 


125. 00000 1.40000 - 00000 - 00000 - 00000 
214 .14E-01 .00 -00 
1414 .12E-01 .00 - 00 
1519 .51E-02 .00 -00 
10 §-.94E-02 .00 00 
12 3-.91E-02 .00 -00 
11.0000 


- 000.000. 000. 000. 000.000. 000. 000.000. 000. 000. 000. 000. 000.000. 000.000. 000.000. 000 
-000.010.010.010.010.010.010.010.010.005. 005.010. 010.010. 010. 000.000. 000. 000. 000 
«000.010.010.010. 010. 010.010. 010. 010. 005. 005. 010.010. 010. 010. 010.010. 010. 000. 000 
- 600.010. 010. 010. 010. 010.010.010.005. 005.010.010.005. 005. 010. 010.010.010.010. 000 
000.000. 000. 000.000. 010.010.010.005. 005. 010.005.005.005. 005.010.010.010. 010.000 
+ 000.000.0000. 000.000.010.010. 010.000. 000. 010.010. 005. 005. 005. 000. 000.010. 000. 000 
- 000.000. 000. 000. 000. 010. 010. 010. 000. 000.020.010.005. 005. 005. 000.000. 010. 000. 000 
- 000.000. 000. 000.000. 010. 010.020. 020. 020.020. 005. 005. 012.012.010.010. 010. 010. 000 
- 000. 060. 060. 020.020. 020. 010.020. 020. 005. 005. 005. 010. 012.012.010.010. 010. 010. 000 
- 000.060. 060. 060. 060.020.005.005. 005.005.010.010. 010. 012.012.025.025. 025.020. 000 
- 000. 060. 060.060. 060. 060. 020.020. 020. 010.010. 010. 025.025. 025. 025.025. 025. 025. 000 
- 000. 060. 060. 060. 060. 060. 060. 060. 060. 020. 010.025. 025.025.025.025. 025.025.025.000 
. 000.060. 060. 060. 060. 060. 060. 060. 050. 020. 020. 025. 025. 025. 025. 025. 025. 025.025. 000 
- 000. 060. 060. 060. 060. 060. 060. 060. 060. 050. 020.025.025.025. 025. 025.025. 025.025.0006 
- 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000.020.020.020. 020.000. 000. 000. 000. 000 
000.000. 000. 000. 000. 000. 000. 000. 000. 000. 000.020.020.020. 020. 000.000. 000.000. 000 
- 000. 000. 000, 000. 000. 000. 000.000. 000. 000.000.015.015. 015. 015. 000. 000.000. 000. 000 
- 000. 000.000. 000.000. 000. 000. 000. 000. 000.000.010.015. 015. 015. 010.006. 006. 006. 000 
- 000. 000.000. 000. 000. 000. 000. 000. 000. 006. 006. 006. 006. 006. 006. 006. 006. 006. 006. 000 
- 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000. 000 
1.10000 

Oo. GO. GO. O. OG. O. O. O. O. GO. O. O. DO. O. GO. O. DO. GO. OD. OL 
O.41.41.41.432.43.45.47.50.50.45.23.12. 5. 2. 0. O. O. O. O. 
0.43.43.43.44.45.437.51.60.72.66.40.20.13. S&S. 5. 2. 2. G O. 
O.44.44.44.46.46.47.52.59.62. 58.42.29. 20.15.10. 7. 7%. 6. O. 

O. O. O O«. 0.42.43. 45.54. 57.55. 36.35.25. 20.16.15.17.17. 0. 

0. O. O. O. 0.40.40.43. 0. 0.50.44.35.29.24. 0. 0.24. O. O. 

0. O. O. O«. 0.44.44.40. 0. 0.47.43.36.32.27. 0. 0.30. GO. O. 

Oo. 0. 0. O. 0.48.45.45.49.51.49.43.38.34.33,33.34.35.36. 0. 

O.51.51.50. 50.50. 47.49.46.47.46.44.40.37.37.37.38.10.40. 0. 

0.51.52. 54.52.52.49.47.47. 46.46. 44.42,42.42,.41.42,43.45. 0. 
0.52.57.57.55.52.50.49.49.47.46.44.44.44.44.45.47.47.50. 0. 

0.53.57.57.55. 54.52.50.49.48.47.47.46.47.47.48.50.52.54. 0. 
0.56.58.58.54.52.49.48.48.49.489.49.47.47.48.49.54.55.57. 0. 
0.59.59.58.54.52.49.49.48.47.46.47.16. 46.18.49.56.57.60. 0. 

Oo. 0. 0. O. O. O. O. O. O. GO O.53.53.54.55. 0. 0. O. 0. O. 

0. Oo. GO. O OO. O O. O O. DO. 0.55.56.60.58. 0. 0. GO. O. O. 

0. 0. O. O. GO. GO. O. GO. 0. GO. O.57.59.60.62. ©. O. O«. GO. O. 

0. O. 0. 0. 0. GO. 0. O O. O O.57.60.63.65.63.60.58.56. 0. 

O. Oo. O. O. GO. O. O. 0. O.54.55.57.60.63.65.63.60.59.56. 0. 

0. O. O. O. 0. GO O. O. 0. O. 0. O. 0. O. 0. O. O. O. DO. OD. 
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APPENDIX I 


SAMPLE OF BIOPLUME II OUTPUT, EDITED 





BIOPLUME il 
CONTAMINANT TBANSPORT UNDBS THE INFLUENCR CF OXYGEN LiMiTBD BLODBGRADATICN 
run-43 Based on 06 May Data, witn $0 Apr - 06 May puaping data 


TSP aU T DATA 
GBID DESCELPTORS 


Ni \NUMBBR OF COLUMNS} 
NY (NUMBER OF ROWS) 
(DBL (i-DISTANCK iN PERT) 
YOBL \Y-DISTANCR IN FBET) 
TIME PARAMETERS 


NTIM «oA MAY. NO. OF TIME STEPS} 
NFMP NO. OF PUMPING PERIODS 


PINT PUMPING PEBIOD IN YEARS} 020 
TIME TIMB INCREMENT MULTIPLIER) 00 
TINIT CLWETIAL THR STEP IN SEC.) 0. 
HYDEOLOGIC AND CHBMICAL PARAKETBRS 
3 (STORAGE COBFFICLENT) = 000000 
POROS  { ERFBCTIVE POROSITY) > 250 
BETA LONGITUDINAL DISPEBSIVITY) = 2.0 
DLTBAT ‘RATIO CP TRANSVERSE TO 
LONGITUDINAL DISPBRSIVITY) = 10 
ANFCTR (RATIO OF T-¥¥ 70 T-2X) = 1.000000 
BXECUTION PABAMETBRS 
NITP (NO. OF ITEBATICN PARAMETERS} = ? 
TL CONVERGENCE CRITERIA - ADIP) = 0010 
“TMA (MANLNO.OP I[TBEATICNS - ADIP] = 100 
_ CRLDIS \MAX.CELL SISTANCE PER HOVE 
QF PARTICLES - 4.0.0.) : 1.000 
WPMAY (MAX. NO. OF PARTICLES) = $000 
NPTPND (NO. PARTICLES FER NUDE) : q 


PROGRAM OPTIONS 


NPNT TIME STBP INTEEVAL FCR 

COMPLETE PSINTOUT} 
NPNTMV (MOVE INTERVAL FOR CHEM. 

CONCENTRATION PRINTOUT) = 0 
NPNTVL (PRINT OPTION-VELOCITY 

Q=NO; L=2PIRST TIME STEP; 

22ALL TIMB STEPS) > 0 
NPNTD (PRINT OPTION-DISP.COEF. 

Q=NO; 1=PIRST TIMB STEP; 

2zALL TIME STEPS) : 0 
NUMOBS (NO. OF OBSBRVATION WBLLS 


— 


FOR HYDROGRAPH PRINTOUT) = 0 
NBEC = (NO. OF PUMPING WELLS} = 5 
NCODES (FOR NODE IDBNT.) > J 


138)8) 





26 


NPNCHY «PUNCH VELOCITIES) > Q 
NPDBLC (PRINT OPT.-CONC. CHANGE} = 0 
REACTION TERMS 
DK (DISTRIBUTION CORFPRICIENT) = .12500B+03 
HOB (BULK DENSITY OF SOLIDS) = .14000B+01 
RF (RETARDATION FACTOB) =, TOL00B+03 
THALE = (HALF LIFR OF DBCAY,IN SBC)=  .00000B+00 
DECAY DRCAY CONSTANT=LN Z/THALF)=  .00000B+00 
DBCAY TERMS 
DBC] !ANARBOBIC DECAY COBFR. | = .000008+00 
DEC2  BBABBATION DECAY CORFF.) = .000008+00 
STEADY-STATE PLOW 
TIMB INTRRVAL (IN SBC} POR SOLUTR-TRANSPORT SIMULATION = .G3L15E+06 


LOCATION OF PUMPING  «2LLS 


1  Y BATE(IN CPS} CONC.  CONC(02) 


ae 0140 00 00 
14 14 0120 00 .00 
18 19 0051 00 00 
i - 0094 00 00 
re ot -,0091 00 00 


ABBA OF CNB CELL = 100.0 
Y-Y¥ SPACING: 

19.000 

10.000 


TRANSMISSIVITY MAP (FTERT/SBC} 


Q.0CE+00 O.008+00 0.002+00 0.008+00 0.008+00 0.00B+00 0.008400 0.008+00 0.008+00 0.002400 
V.008+00 0.008400 0.008400 0.00K+00 0.00E+00 0.008+00 0.002400 0.008+00 0.00B+00 0.002400 
Q,002+00 1.00B-02 1,00B8-02 1.008-02 1.00B-02 1.00B-02 1.00E-02 1.002-02 1.008-02 5.008-03 
5.00B-03 1.00B-02 1.008-02 1.00B-02 1.00B-02 0.00B+00 O.00E+00 0.008+00 0.00F+00 0.008+00 
0.002+00 1.008-02 1.G0B-02 1.00B-02 1.00B-02 1.00B-02 1.00B-02 1.008-02 1.008-02 5.008-03 
5,0CB-03 1.00B-02 1.,008-02 1.008-02 1.008-02 1.00B-02 1,608-02 1.008-02 0,002+00 0,008+00 
Q.00E+00 1.008-02 1.002-02 1.00B-02 1.00B-02 1.008-02 1.002-02 1.002-02 5.008-03 5.008-03 
1.00B-02 1.008-02 5.008-03 5.002-03 1.00B-02 1.00B-02 1,008-02 1.00B-02 1,008-02 0.008+00 
0.00B+00 0.002400 0.002400 0.008400 0.00E+00 1.008-02 1.002-02 1.002-02 5.002-03 5.008-03 
1,00B-02 §.002-03 5.00R-03 §.008-03 §.00B-03 1.008-02 1,008-02 1.008-02 1.008-02 0.002400 
0.00E+00 0.00B+00 0.002400 0.002+00 O0.00E+00 1.00E-02 1.00E-02 1.008-02 0.00B+00 0.008+00 
1.008-02 1.008-02 5.008-03 5.008-03 5.00E-03 0.008+00 0.002400 1.008-02 0.00B+00 0.008+00 
0.008400 0.00B+00 O0.008+00 0.008400 0.008+00 1.008-02 1.008-02 1.008-02 0.00B+00 0. 00B+00 
2.00B-02 1.00B-02 5.00B-03 5.002-02 5.00B-03 0.00B+00 0.002+00 1.008-02 0.00B+00 0.008400 
0.008+00 0.00R+00 0.008+00 0.00E+00 6.008+00 1.008-02 1,008-02 2.008-02 2.00B-02 2.00B-02 
2,00B-02 §.008-03 5.008-03 1.208-02 1.208-02 1.008-02 1,008-02 1.008-02 1.008-02 0,008+00 
0.008+00 6.008-02 6.00B-02 2,00B-02 2.008-02 2.00B-02 1,00B-02 2.002-02 2.00B-02 5,008-03 
5.00R-03 §.008-03 1.008-02 1,208-02 1.208-02 1.008-02 1.,00B-02 1.008-02 1.00B-02 0,008+00 


baja) i! 





),.008+00 
1.008-02 
5, 008+00 
1008-02 
3, 00B+00 
1008-62 
7, O08+90 
2. 00B-02 
6, 00B+00 
2,00B-02 
0.008+00 
).00B+00 
0.00B+00 
),00B+00 
0.008+00 
0.00B+00 
).008+00 
J. 008+00 
0.00B+00 
6 .00B-03 
), 00B+00 
0. 00B+06 


008 +00 
5 0B-02 
 00B+00 
008-02 
-008+00 
008-03 
) .008+00 
0, 002+00 


no KK CO r—M Oo moescCt 


ANUIFER THICKNESS (FT) 


Alm .0 


ned. 


ne 5.1 


nS. 


Sil 


Sat 


Hod 


00 9.3 §.7 


ah 


5.9 


5.8 


5.8 


208-02 
008-02 
00B-02 
50R-02 
008-02 
50R-02 
008-02 
50B-02 
 0OR-02 
50B-02 
.008+00 
-00R-02 
00B+00 
COB-62 
-O0B+00 
5OB-02 
 00B+00 
508-02 
O0B+00 
00B-03 
008+00 
008400 


5.0 


502 


5.8 


5.4 


5.4 


5.2 


tee 


0OR-02 
COR-02 
VOR-02 
508-02 
O0E-02 
50R-02 
008-02 
50R-02 
0OR-02 
50R-02 
 OOB+00 
 VOR-02 
 00B+00 
 0OR-02 
 00R+00 
5OR-02 
OOB+00 
50E-02 
 00B+00 
008-03 
008+00 
OOF +00 


4.5 
4.7 
4.7 
4.3 
4.0 
4.4 
4.5 
4.1 
4.3 
5.0 
5.2 
43 


4,9 


6.00B-02 
1.208-02 
5. 0UR-02 
2.508-02 
6.00B-02 
2.50B8-02 
€.008-02 
2.50B-02 
6.00R-02 
2508-02 
0.00B+00 
2.008-02 
0.00R+00 
2. 005-02 
0.008+00 
1.50E-02 
).008+00 
1.508-02 
0.008+00 
§.00B-03 
0). 008+00 
0.008400 


Fs 
Suc vom 


Mee 


om 
a 


Lee iee eh, 


2.00R-02 
2. 508-02 
€.00B-02 
2.50B-02 
6.00R-02 
2. 5CR-02 
6.008-02 
2. 50B-02 
6. 00R-02 
2. 508-02 
0.002+00 
0.008+00 
0 .00R+00 
0.008+00 
0.008+00 
0.008+00 
0. 00B+00 
1. 00B-02 
0.008+00 
6.00B-03 
,002+00 
0.002+00 


Shall es 
Gers gl) 
6.2 mee 


Ded were 


Se 
4.7 4.6 
1.6 496 
4,7 eb 
4.80 4,7 


4.9 4.8 


cro on 


cro om © oO Fo Fe 


6. 
a 
0. 
0. 
0. 
0. 
0, 
0. 
0. 
6. 
0. 
6. 
0. 
0, 


QUE-03 
308-02 
QOB-02 


J0B-02 
-O0B-02 
.50B-02 


00E-02 
50B-02 
OOR-02 
50B-02 
)0B+00 
008+00 
00E+00 
00B+00 
002+00 
0OR+00 
008+00 
00-03 
COER+00 
00R-03 
00B+00 
00E+00 


co 
= 


cc 
ae 


ars 


6) 


3.6 


3.8 


4.0 


4.2 


4,4 


4,6 


4.1 


4.6 


UUR-03 
UR-02 
OOR-02 
50B-02 
008-02 
50R-070 
QOR-02 
502-02 
VOB-02 
50R-02 
008 +00 
OOE+00 
OUE+00 
 O0R+00 
 OOB+00 
008+00 
00B+00 
008-03 
 008+00 
002-03 
0.00B+00 
0.00B+00 


Oo © OD es tn 


ronmrmreraoece coe re Nh oH fs 


0. U0E+00 
0. 00B+00 
). VOB +00 
0.00B+00 
0. U0B+00 
0, UUB+00 
0. OUR+G0 
§ UUB-09 
0, U0B+00 
b, OUB-03 
0. 00B+00 
0. SOB+00 


ec 
o> 


3.7 


4,1 


4,5 


4.8 


4,9 


5.0 


5.4 


5.6 


nay] 


& U0B-03 
U.0OE+00 
1, 0uB-02 
0, 00B+00 
2.00B-02 
0), OUB+00 
U0B-02 
O0B+00 
008-62 


VOR +00 
OOB+00 
0. 0OB+00 
0, 00B+00 


a co cet Co COC 


ane 


5.8 


5.7 


0 
) QOB+00 
0 
0 


), OUR +00 
OUR +00 
OOE+00 
QUE +00 
 OUE-03 
QOR+00 
008400 
 QOB+U0 


3.6 


4.5 


5.0 


5.4 


3.7 


6.0 





0 
iO 

0 
0 

0 
0 

0 
0 

0 
0 

0 
0 


all lt 
? .0 
0 .0 
0 ad 
0 9 
a) 


Dd 0 
0d 
0 .0 
0 wl 
0 0 
0 0 


DIFFUSE BBCHARGE AND DISCHARGE (FT/SEC) 


0.008+00 
0. 008+00 
0. 008+00 
2. 108-09 
0.00B+00 
2. 10B-09 
0.008+00 
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0.008+00 
9.00B+00 
0. 00B+00 
0. 00B+00 
0.00B+00 
0. 00B+00 
0.00B+00 
0.008+00 
.00B+00 


1.00B+00 


Uae 
0. 
a. 
VE 
0, 
0. 0. 
0 0 
t. 0. 
0. 0. 
0, 0. 
0. 0. 
do 0 
0 0 


0. 008+00 
0.008+00 
0.008+00 
0.008+00 
0.00E+00 
0.008+00 
0. 008+00 
0.008+00 
0.008+00 
0.008+00 
0.90B+00 
0.008+00 
0.008+00 
0.008+00 
0.008+00 
0.008+00 
0. 008+00 
0.008+00 
0.008+00 
0.00E+00 
0. 008+00 
1.008+00 
0.008+00 
0. 008+00 


0,008+00 
0.002+00 
0, 008+00 
0.008+00 
0.00R+00 
0.008+00 
0.008+00 
0.008+00 
0. 00B+00 
0.008+00 
0.00B+00 
0.00B+00 
0.002+00 
0.008+00 
0.008+00 
0.008+00 
0.00R+00 
0. 008+00 
.008+00 
0.008+00 
0.002+00 
1.008+00 
0.008+00 
0.00B+00 


0.00B+00 
0, 00B+00 
0, 00E+00 
0.90B+00 
0.00R+06 
G.00E+00 
0.008+00 
0,00B+00 
0, 008+00 
0. 008+00 
0. 00B+00 
).00B+00 
0.00B+00 
0.008+00 
0. 008+00 
0.00R+00 
0.00E+00 
0.00B+00 
0.008+00 
0.00B+00 
0.008 +00 


1.00B+00 


0.00B+00 
0, 00B+00 


}, 00B+00 
0. 008+00 
‘), 0OR+00 
0,00R+00 
§008+00 
), O0B+00 
0, 00B+00 
}, 00B+00 
0. 00B+00 
0, O0R+00 
), 008+00 
0. 00R+00 
0.00B+00 
0.008+00 
0. 00B+00 
), 0UR+00 
), 008+00 
0.008+00 
0. 00B+00 
0. 00B+00 
0, 008+00 
1. 00B+00 
0.005+00 
0, 008+00 


0. 00B+00 
0. 00B+00 
0. 00B+00 
0. OUE+00 
), OOR+00 
).008+00 
0. 00B+00 
0.00K+00 
0. VOB +00 
0. UUE+O0 
0, 0OB +00 
0.00R+00 
0. 00E+00 
0.00R+00 
0.008 +00 
0. 00E +00 
0.00B+00 
0. 00R+00 
0), UOB+00 
0), 00R+00 
0.008 +00 
1. QOB+00 
). DUB +00 
0. 00B+00 
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0, 008 +00 
0. 008+00 
0. 008+00 
).00B+00 
0. 0OB+00 
0, Q0E+00 
0. U0R+00 
0.008 +00 
0. 00R+00 
0. 00R+00 
0, 00E+00 
0.008 +00 
0, 00B+00 
0.00B+00 
0. 00R+00 
0.008400 
0. 00B+00 
0.002400 
0.008+00 
0). O0R+00 
1.00B+00 
0.008400 
0. 00B+00 
0.008+00 





0. 
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0, 
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0. 
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0. 
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th 
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0. 
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0. 
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HBAD DISTBIBUTION - &GW 


NUMSBB OF TIME STEPS 
TIME (SECONDS) 
TIME( DAYS) 
TIMB( YEARS) 


0 = 0000000 
0000000 

0 0000000 
0000000 

0 0000000 
- 6000000 

0  .0000000 
0000000 

0 «=. 0000000 


§.2100000 = 3. 


0 §.1300000 
0000000 

0 «=. 0000000 
- 0000000 

0 0000000 
- 0000000 

0 =. 0000000 
0000000 

0  .0000000 
- 0000000 

2 0000000 
- 0000000 

0 «=. 0000000 
- 0000000 

0 «=. 0000000 
- 0000000 

0 = §=—_- 0000000 
0000000 

0  .0000000 
0000000 


0000000 


-0000000 


0000000 


0000000 


0000000 


-0000000 


0000000 


0000000 


73000000 
1300000 
3.1100000 


. 0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


- 0000000 


- 0000000 


0000000 


0000000 


. 0000000 


-0000000 


-0000000 


0000000 


-0000000 


0000000 


- 0000000 


0 
00000 
- 000008+00 
-000008+00 
0000000 
0000000 
. 0000000 


. 0000000 


= 


3700000 


aD 


0500000 


0000000 


0000000 


0000000 


. 0000000 


0000000 


0000000 


0000000 


0000000 


. 0000000 


— 


—4 


0000000 


0000000 


0000000 


0000000 


.5000000 


. 9500000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


~ 


~~ 


te 
i 


0000000 


0000000 


. 0000000 


0000000 


. 1000000 


- 8500000 


0000000 


(000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


7, 9000000 


. 0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


0000000 


a 


0000000 
0000000 
0000000 
0000000 
. 1000000 
0000000 
0000000 
. 0000000 
0000000 
. 0000000 
0000000 
0000000 
0000000 
0000000 


-0000000 


oD 
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0000000 
0000000 
0000000 
. 0000000 
2900000 
0000000 
9000000 
. 0000000 
. 0000000 
0000000 
- 0000000 
0000000 
0000000 
. 9000000 


- 0000000 





ee 
tl lily bbe 


Be 


0 0000000 
0000000 

9  ~—.0000000 
0000000 

0 0000000 
0000000 

0 0000000 
0000000 

0 0000000 
0000000 

d 0000000 
. 0000000 

0000000 
0000000 

0 0000000 
0000000 

0 0000000 
0000000 

: 0000000 
. 0000000 

0 0000000 
. 0000000 

0 0000000 
0000000 

0 0000000 
. 0000000 

0 .0000000 
0000000 

0 0000000 
0000000 

0 .0000600 
0000000 

0 . 0000000 
0000000 

0 0000000 
0000000 

0 0000000 
. 0000000 

0 0000000 
0000000 

0 0000000 
0000000 

0 0000000 
0000000 

0 2.9300000 

2.0300000 

0 0000000 
0000000 

0 =. 0000000 
0000000 


- 0000000 
0000000 

0000000 
0000000 

0000000 
4000000 

0000000 
0000000 

-9000000 
0000000 

0000000 
0000000 

0000000 
0000000 

. 0000000 
. 0000000 

0000000 
0000000 

0000000 
. 0000000 

. 0000000 
- 0000000 

0000000 
0000000 

-0000000 
0000000 

- 0000000 
-0000000 

0000000 
. 0000000 

0000000 
- 9000000 

0000000 
0000000 

- 0000000 
0000000 

. 000000 
0000000 

0000000 
, 0000000 

. 9000000 
0000000 

0000000 


2.0300000 


2.0300000 
0000000 

0000000 
0000000 

-0000000 
0000000 


ILTEBATION PARAMETERS 


2597262-02 


. 0000000 
- 0000000 
. 9000000 
. 0000000 
. 0000000 
0000000 
0000000 
-0000000 
0000000 
- 0000000 
0000000 
. 0000000 
0000000 
0000000 
0000000 
. 0000000 
0000000 
0000000 
0000000 
- 0000000 
. 0000000 
. 0000000 
2.0300000 
0000000 


0000000 


. 9000000 
0000000 
0000000 
0000000 
. 0000000 
0000000 
. 0000000 
0000000 
. 0000000 
. 0000000 
0000000 
0000000 
. 0000000 
0000000 
. 0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
2.0300000 
0000000 


0000000 


ss 


0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
. 0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
. 0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0300000 
0000000 


0000000 


i 7 


1000000 
0000000 
. 0000000 
. 0000000 
9000000 
- 0000000 
. 0000000 
4000000 
0000000 
- 0000000 
0000000 
. 0000000 
0000000 
0000000 
9000000 
0000000 
0000000 
. 9000000 
0000000 
0000000 
0000000 
0000000 
0300000 
- 0000000 


0000000 


cs 


. 0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0000000 
0300000 
0000000 


0000000 


cr 


163 


9000000 
. 0000000 
0000000 
0000000 
- 1000000 
0000000 
0000000 
. 9000000 
- 9000000 
0000000 
1000000 
0000000 
9000000 
0000000 
. 9000000 
0000000 
0000000 
0000000 
- 0000000 
0000000 
0000000 
0000000 
0300000 
0000000 


0000000 
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-1005358-02 
1889492-01 
9096348-01 
137459 
370754 
1.00000 
000000 
000000 
000000 
000000 
. 000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 


CONCENTRATION OF CONTAMINANT 


NUMBBR OF TIMB STRPS = 0 
TIME(SBCONDS) = 00000 
CHEM. TIMB(SECONDS) =  .000008+00 
CHEM. TIMB( DAYS) = —- 000008+00 
TIMB(YEARS) == =. 000008+00 
CHEM. TIMB( YEARS) = .000008+00 
NO, MOVES CCMPLETED = 0 


ioe gf 0 0 Odo 9 0 ete tS 0 eee 80 re | 0 
» 9 et! ® 0 Oo @O 8 © 100 100 ee era v0 0 
0 0 o 90 0 Oo 0 0 0 0 100 100 300 oO 90 0 0 0 0 
> 2 0 0 0 100 100 100 100 500 0 O 0 C0 
7 0 Oo O 0 0 © 100 100 100 100 100 50 OO 0 0 ’ 0 0 
0 0 0 O 0 8000 5000 100 0 0 100 200 1000 en | o 9 0 
0 Oo oO 9 010000 6000 4000 §=6©0 «©6200» 300 1000 = 100 art: 0sr out 
0 0 0 0 010000 6000 4000 100 100 200 300 2000 100 0 0 ® od 90 
0 0 0 §00 150010000 6000 4000 500 200 200 1000 3000 200 100 oD » 90 0 
0 9 §00 2000 25001000010000 5000 5000 5000 5000 5000 5000 2000 1000 100 Q ) 


0 1000 2000 70001000015000150001500015000100001000010000 7000 5000 1000 1000 en 


@Peeorocascnocioeoceororee ere em eo a cr ec ec oe fc ef & & 





G —-& 3000 6000 75003000030000300003000030000200001500010000 


J 0 -2500 75001000030000300003 00003000030000300002900010000 


0 
0 
) 90 
d 60 
0 
0 9 
0 
i ¢ 
0 0 
0 
ie 0 
U 


0 


0 


0 


) 
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CONCENTRATION OF OXYGEN 


NUMBER OF TIME STBPS 


TIME( SECONDS) 
CHEM. TIMR( SECONDS) 
CHBM. TIMB (DAYS) 

TIMB( YEARS) 
CHEM. TIME ( YEARS) 


NO. MOVBS COMPLETED 


0 


0 


gooodoncnrcqocnroocoocooocococ ew oc eco OO lS 
= 


0 


0 


0 


0 


0 


0 


) 


0 


d) J 
0 9 
0 0 
@ 0 
0 0 
0 0 
00000 
-000008+00 
-00000E+00 
-000008+00 
-000008+00 
0 
0 9 
0 0 
0 9 
0 0 
0 oO 
0 9 
0 0 
0 0 
0 2 
9 90 
0 0 
0 9 


0 


0 


0 


Y 806010000200003000030000300003000030000300002000012000 


0 1000 
0 300 
0 300 
0 200 
0 100 
0 9 
0 0 
0 0 
0 0 
0 40 
0 9 
0 9 
0 9 
0 0 
0 90 
0 0 
0 9 
0 0 


7000 Sy00 
$000 800 
8000 7000 
2000 6000 
300 © 300 
300 300 
200 © 300 
200 300 
0 od 

y 

i) Q 

0 i 

0 0 

0 9 

0 0 

0 0 

? 9 

b 9d 

0 ) 

0 0 


Lyd 


) 


v 


3000 2000 
$500 Ju0d 1000 
6000 4500 2000 

10 U 

300 y 

300 

200 ) 

100) 

en) 

i) y 

oO VW 

0 0 

do 0 

iad 

0 90 

0 9 

db 0 

0 2 

0 0 

0 

do 





0 
0 ) v ) ) v 
0 
0 0 0 v ) u 0 
0 
0 i] 1) i) y) ) g ) 
0) 
y 0 ) 0 0 i} 0 0 
0 
0] ) 0 1) i} 0 j ) 
‘ 
Q Q 0 0 0 i} 
0 
d 1} ) i) 0 } 0 
) 
0 0 0 0 ) 0 0 i] 
i) 
0 
N= l 
NUMBER OF [TSRATIGNS = 35 
LHEAD DISTBIBUTION - &CW 
NUMBB2 OF TIME STEPS = l 
TIMB(SBCONDS) = .G3115R+06 
TIME (DAYS) =, T3050E+01 
TIMB(YEABS} = .200008-01 
i] 0000000 9000000 0000000 
0000000 9000000 
0 0000000 0000000 (000000 
0000000 0000000 
q 0000000 =. 7.3831925 = 7. 4263873 
§.1809111 &.1460156 
0 6.1155542 8.0819247  3.0314900 
9000000 0000000 
) 0000000 =7.3000188  7.3700156 
%.2099778 §=8. 1299917 
0 §8.1299739 3.1099609 8.0499727 
0000000 0000000 
Q 0000000 «=7.3547790 = 7. 4095387 
7,8802499 7.9422036 
0 7.7698848 7.7932957 17.7276529 
6.9480738 9000000 
0 0000000 0000000 0000000 
7.5896384 7.6128273 
0 7.4582205 7.4190843 7.3194518 
6.8789457 0000000 
0 0000000 0000000 0000000 
0000000 0000000 
0 7.0695949 7,0804853  6.9235490 
0000000 0000000 
0 0000000 0000000 0000000 


0000000 


0000000 


0000000 
0000000 
7,5459748 
7.9625726 
7.50001 20 
7.9499792 
7,5038195 
7.6185467 
0000000 
1.207471 
0000000 
6.7263366 


0000000 


1 
: 


1 


1 


1 


1 


6 


0000000 
0000000 
T115460 
9062485 
6999913 
8499455 
6018982 
5136005 
0000000 
1655110 
0000000 
6271500 


. 0000000 


0000000 


0000000 


— 


»§886930 


0000000 


a-84 


8999691 


~ 


«4849065 


—> 


6018822 


— 


4394826 


—= 


2141964 


—3 


. 1898508 


6.7736825 


0000000 


62986533 


2000000 
0000000 
§.0545849 
2000000 
3.0999536 
7.2651767 
7.6914830 
7.1695662 
1.26705 47 
70563626 
§. 8082244 
0000000 


63501179 
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0000000 


0000000 
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EN esh| 


0000000 


ao 


2899122 


= 
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= 


1970542 


= 
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= 


9543290 
§.8096729 
68420441 
6. 2860237 


§.3635868 
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§.7699750 
. 0000000 
9900000 
0860698 = 6 
6.4367204 
0888331 
. 0000000 


= 


58379798 
§391588 
0000000 


5.244314) 
8040454 
. 000000 
cee) an 
4.9192468 
7203708 
0000000 


Ce el oe el — cD 


46918913 
6498051 
0000000 


as 


4.5580797 
5939112 
0000000 


4.4921107 
5623282 
. 0000000 
0000000 
0 0000000 
0000000 
0 0000000 
. 0000000 
0 —. ¥000000 
. 0000000 
0 0000000 
0000000 
0 0000000 
0000000 
0 0000000 
0000000 
0 0000000 
© .0553932 
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0 = 0000000 
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17045844, 


toeelgl =. 


6.9702717 


- 0000000 
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0000000 


48031795 
9147063 
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5.1177468 
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. 0000000 
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0000000 
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0000000 
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0000000 
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0000000 


1000000 


. 0000000 
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0000000 


. 0000000 


0000000 
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0000000 
9300000 
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0000000 


. 0000000 


0000000 
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0000000 
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48053841 
4,9963909 
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4,7079305 
4.7960449 
45231246 
4.637045] 
4,3202082 
4.638991 
40758496 

0000000 
2 8043937 

0000000 
3. 4969877 

0000000 
2,1369598 

0000000 
2.7181048 

. 0000000 
2.0300590 


0000000 
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on 
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_ 
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en 
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ay 
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=~ 
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—_ 
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~_ 
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oo 


4970209 


a 
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en 
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0000000 


4662201 


en 


. 1000000 
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0000000 


2.683213 
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= 
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— 
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— 


13345599 
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. 0000000 


0000000 
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0000000 
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tc 


3194397 
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0000000 


1000000 
59299054 
3.1333073 
55544027 
49636378 
51557264 
48003386 
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4$,7039051 
48261190 
4,6044355 
47965447 
45177961 
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4,4609804 
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0000000 
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- 0000000 
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0000000 
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0000000 
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. 9000000 
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0 CUMULATIVE MASS BALANCE -- (IN PT333) 


0 


RECHARGE AND INJECTION 


PUMPAGE AND B-T WITHDRAWAL 
NBT PUMPAGE 
WATER RELEASE FROM STORAGE 
INTO AQUIFER 


CUMULATIVE 


LEAEAGE 


LEAKAGE OUT OF AQUIFER 


CUMULATEVE 


NET LEARAGE 


MASS BALANCE RESIDUAL 


ERROR 


(AS PBSCRNT) 


- 116762405 
136468405 
T9GG2B+O4 
000008 +00 
CObLIBHOS 

- 156312405 
197 98B+04 
10.541 
29876E-01 


Q BATE MASS BALANCE -- (IN C.P.S., 


0 


0 


LEABAGE INTO AQUIFER 
LEAKAGE OUT OF AQUIFER 


NET 


LRAEAGE 


(QNET) 


RECHARGE AND INJECTION 


PUMPAGE AND B-T WITHDRAWAL 


NET 


WITHDRAWAL 


(TPUM) 


STHO9E-01 
- 24T66E-01 
126438-01 
-.18500B-01 
»311268-01 
= 126262-01 


STABILITY CRITERIA --- 4.0.¢. 


FLUID VELOCITIES 
YHAY = 3, 448-03 VAY = 1.248-03 
VMEBD= 2, 518-03 VMYBD= 1.1 78-03 
EFFRCTIVE SOLUTE VELOCITIES 
VAT = 4,908-06 VMAY = 1.778-06 
VMEBD= 3. 588-06 VMYBD= 1. 678-06 
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) EFPBCTIVE CXYGEN VELOCITIES 
WAT 3 2. 448-03 VMAY = 1.248-03 
WMXBD= 2.51B-03 VMYBD= 1.178-03 
U THY (MAY. INJ.) = 3669.0 
SEMV (CELDIS) = 2909.0 
Q TIMV = 2.918403 NTIMV = = 216 NMOV = 217 
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NO. OF PARTICLE MOVES REQUIRED TO COMPLETE THIS TIMB STEP = 217 


CONCENTRATION OF CONTAMINANT 
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NUMBEE OF TIME STEPS 


DELTA T = 631158406 

TIME(SBCONDS) = .63115B+06 

CHEM. TIMB(SECONDS) = .S31152+06 

ABM. TIME (DAYS) =. 730508+01 

TIMB(YEABS) = .cOQO0B-01 

CHBM.TIME(YBARS: = .200008-01 
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CHEMICAL MASS BALANCE 
MASS IN BOUNDARIES =. 626178+06 
MASS OUT SOUNDABIES = -.013658+07 
MASS PUMPED IN = ,000002+00 
MASS PUMPED OUT = =, 124178409 
MASS LOST ¥. BIODEG. = .285648+05 
MASS LOST BY BADIO. DCY=  .0000CE+00 
MASS LOST BY ANABB. DCY= .000008+00 
MASS LOST BY BEABR. DCY= .000002+00 
MASS ADSOBBBD ON SCLIDS= .914832+411 
INITIAL MASS ADSORBED = .91455E+11 
INPLO® MINUS OUTFLOW = -.125682+09 
INITIAL MASS DISSOLVED = .{3065B+09 
PRESENT MASS DISSOLVED = .130622+09 
CHANGE MASS DISSOLVED = -.359208+04 
CHANGER TOTL.MASS STORBD= -.251032+07 
COMPARE RESIDUAL @ITH NBT FLUX AND MASS ACCUMULATION: 
MASS BALANCE BRSIDUAL = -.ic3162+09 
ERROR (&3 PERCENT) = -.193608+05 
COMPARR INITIAL MASS STORED WITH CHANGE IN MASS STOBBD: 
ERROR (AS PERCENT) =. 134308+00 
CONCENTRATION OF OXYGEN 
NUMBER OF TIME STBPS = 1 
DELTA T = ASL152 +06 
TIMB(SECONDS) = .631158+06 
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CHEMICAL MASS BALANCE FOR OZYGEN 


MASS IN BOUDARIBS = 942498405 
MASS OUT BOUNDABIRS = -,189368+05 
MASS PUMPED IN = —. 000008+00 
MASS PUMPED OUT = -, 127448406 
MASS LOST W. BLIODBG. = .856938+05 
INFLOW MINUS OUTFLOW = -.112128+06 
INITIAL MASS DISSOLVED = .00000B+00 
PRRSBNT MASS DISSOLVED = .97627B+04 
CHANGE MASS DISSOLVED = .95455B+05 


CHANGE TOTL.MASS STORBD=  .95455B+05 

COMPARB RESIDUAL #ITH NET FLUX AND MASS ACCUMULATION FOR OXYGEN: 
MASS BALANCE BESIDUAL = -.20758B+06 

BRROR (AS PRRCBNT) = -,606098+03 
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